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ABSTRACT 


This  report  describes  the  statistical  performance  of  several  radar-based  adap¬ 
tive  detection  schemes  in  both  stationary  and  nonstationary  noise  and  interference 
environments.  The  detectors  under  study  must  be  able  to  correctly  determine  the 
presence  of  a  target  in  a  range  gate  with  a  high  degree  of  probability  given  that  the 
probability  of  raisclassification  is  a  fixed  small  value.  The  hostile  noise  environment 
is  assumed  to  consist  of  possibly  time- varying,  spatially  correlated  interference  along 
with  Gaussian  background  noise.  In  a  typical  radar  environment,  the  mean  value  of 
the  returned  radar  signal  and  the  noise  covariance  matrix  are  unknown  parameters; 
therefore,  generalized  likelihood  ratio  test  procedures  were  used  to  develop  deci¬ 
sion  rules  that  meet  the  Neyman-Pearson  criterion.  Three  major  cases  of  interest 
were  examined.  First,  the  single-pulse  test  developed  by  Kelly’  is  reviewed.  The 
multiple-pulse  return  test  case  is  extremely  complicated  and  was  divided  into  dis¬ 
tinct  detector  forms;  noncoherent  and  coherent.  The  performance  of  each  detector 
is  a  function  of  the  signal-to-nolse  ratio,  the  number  of  radar  pulse  returns  used  in 
the  decision  rule,  and  the  quality  of  the  covariance  estimate.  Ao^ea'jion  for 
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1.  INTRODUCTION 


This  section  establishes  the  basic  premises  of  this  investigation.  The  intent  of  this  study  is  to 
quantify  the  statistical  performance  of  several  detectors  derived  from  a  generalized  likelihood  ratio 
(GLR)  test.  A  general  overview  of  the  detection  process  is  presented,  and  the  terminology  and 
assumptions  used  in  the  report  are  defined.  Fbur  separate  cases  are  considered,  and  the  general 
problem  definition,  solution  techniques,  and  associated  results  are  briefly  discussed  for  each. 

1.1  Gleneral  Overview 

This  report  addresses  the  development  and  analysis  of  radar-based  adaptive  decision  rules. 
The  particular  concern  is  to  determine  whether  the  returned  radar  signal  indicates  the  presence 
of  a  target  In  the  range  gate  currently  under  observation.  If  detector  output  exceeds  a  threshold, 
certain  actions  such  as  tracking  can  be  initiated.  Because  the  returned  signal  includes  a  possible  echo 
amplitude  proportional  to  target  radar  cross-sectional  area,  directional  interference,  and  Gaussian 
sky  noise,  a  statistical  approach  to  decision  making  is  taken  to  minimize  the  probability  of  an 
incorrect  decision. 

It  is  important  to  note  that  the  term  “adaptive”  used  in  the  current  context  refers  to  the 
time  dependency  of  the  actual  parameter  values  used  in  the  decision  rule.  Adaptive  does  not  refer 
to  the  time  dependency  of  parameter  value  updates  that  occur  in  algorithms  such  as  the  least  mean 
square  and  recursive  least  squares  (1).  To  expand  on  this  terminology  further,  the  scalar  output  of 
an  adaptive  array  processor  is  defined  as 

yi=w«x,  ,  (I) 

where  Wj  is  the  weight  vector  and  Xj  is  the  data  vector. 

At  each  sample  instant  the  adaptive  algorithms  update  the  weight  vector,  the  value  of  which 
is  determined  according  to  a  specified  cost  function.  This  study  deals  purely  with  raw  input  data 
and  does  not  form  the  adaptive  array  output.  Estimates  of  the  parameters  found  in  the  subsequent 
decision  rules  are  time  varying.  Also,  the  decision  rule  contains  forms  similar  to  the  optimum 
weight  vector  for  the  known  covariance  matrix  case. 

In  almost  all  cases,  this  report  is  occupied  with  the  presence  of  noise  and  interference  that 
causes  errors  in  the  receiver  output.  Fbr  a  system  in  which  the  signal  is  corrupted  by  noise,  it  is 
most  desirable  that  the  receiver  minimizes  the  probability  of  error.  Statistical  decision  theory  can 
be  ised  to  obtain  expressions  for  the  probability  of  error,  l.e.,  the  probability  of  misclasslfying  the 
received  binary  signal,  and  the  probability  of  detection  (Pp)  [2].  The  misclassification  error  that  is 
of  most  interest  to  radar  system  analysts  is  the  probability  of  false  alarm  (Pfa)-  This  error  occurs 
when  the  decision  rule  indicates  that  a  target  Is  present  when  actually  it  is  not.  In  practice,  radar 
systems  typically  require  a  Pfa  on  the  order  of  10“®. 
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Several  additional  levels  of  complication  c<’.i  b-  .id-bed  to  this  scenario.  Assuming  white 
Gaussian  receiver  noise  and  a  multiple-elemen'  (.<:■  aray,  the  probability  density  function 
(pdf)  of  the  received  data  is  uniquely  specifieo  b,  'wo  (.arameters — the  mean  value  vector  and 
the  covariance  matrix.  The  mean  and  covariance  matrix  are  also  referred  to  as  the  “first-  and 
“second-order”  noise  statistics.  If  one  both  parameters  are  unknown  and/or  time  varying,  the 
analysis  can  be  become  ext remel-  diilicJ  if  not  intractable.  This  report  focuses  on  examining 
conditions  where  the  noise  enviroi..'.  ''f  is  nonstationary,  i.e.,  where  the  first-  and/or  second-order 
noise  statistics  are  time  varying.  In  addition,  the  true  values  of  the  parameters  are  unknown. 
Section  1.2  describes  the  general  radar  system  model  used  in  this  investigation. 

1.2  Radar  System  Model 

A  comprehensive  review  of  basic  radar  principles  and  operations  can  be  found  in  sources 
such  as  the  Radar  Handbook  [3].  A  radar  system  scans  a  three-dimensional  region  of  space,  the 
coordinates  of  which  are  azimuth,  elevation,  and  range.  The  search  volume  is  divided  into  segments 
called  range  bins  or  range  gates.  To  detect  a  target,  a  sequence  of  modulated  pulses  is  emitted  and 
the  system  then  listens  for  the  returned  echo.  The  time  between  transmission  and  echo  reception 
corresponds  to  the  total  propagation  delay  that  is  a  function  of  distance.  The  received  data  are 
processed  and  input  to  a  decision  rule  to  determine  whether  a  target  is  present  in  the  range  bin 
under  observation.  If  the  output  of  the  decision  rule  is  greater  than  some  predefined  threshold  with 
a  value  that  is  proportional  to  the  desired  false  alarm  rate,  a  “hit”  is  declared.  In  general,  if  a 
sufficient  number  of  hits  occurs  in  a  range  gate,  further  action  is  taken. 

1.2.1  Array  and  Receiver  Structure 

The  radar  system  under  examination  is  an  array  with  N  elements  that  are  connected  to 
separate  receiver  channels.  The  array  configuration  is  not  specified  but  a  generalized  steering 
vector  is  denoted  as  the  TV-vector  p-  The  steering  vector  defines  the  directional  orientation  of 
a  possible  target  and  is  determined  from  target  position  and  array  geometry.  By  not  assuming 
any  antenna  characteristics  or  special  array  configurations,  the  problem  can  be  analyzed  in  its 
most  general  form;  however,  for  subsequent  simulation  studies  a  particular  array  geometry  will 
be  designated.  The  receiver  section  of  the  radar  system  is  assumed  ideally  to  demodulate,  match 
filter,  and  sample  the  incoming  planar  narrowband  waveform.  The  current  model  is  idealized  to 
concentrate  solely  on  first-order  eflfects. 

If  the  target  is  present  in  the  range  bin  currently  under  observation  and  designated  as  hypoth¬ 
esis  Hi,  the  received  signal  consists  of  Nf,  modulated  pulses  plus  noise.  Under  the  null  hypothesis. 
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designated  as  Ho,  the  received  signal  consists  solely  of  noise.  A  sample  of  the  matched  filter  out¬ 
puts  at  time  instant  f,  j  =  1, . . .  ,Np,  is  denoted  as  the  complex^  A’-vector  Zi-  The  vector  Zi  is 
referred  to  as  a  primary  data  vector.  The  samples  of  the  matched  filter  outputs  are  assumed  to 
be  uncorrelated  with  respect  to  time,  but  there  can  be  time- varying  spatial  correlations  among  the 
multiple  receiver  channels. 

In  radar  applications  the  matched  filter  output,  which  directly  corresponds  to  the  complex 
returned  signal  amplitude,  is  not  known  a  priori  and  as  such  must  be  treated  as  an  unknown  pa¬ 
rameter.  The  analysis  is  restricted  to  the  case  where  the  returned  signal  amplitude  is  a  nonrandom 
variable,  i.e.,  the  target  is  nonfluctuating;  therefore,  Swerling  fluctuating  target  distribution  models 
[4]  are  not  included  in  this  analysis.  The  complex  returned  signal  amplitude  is  assumed  constant 
over  the  Np  observation  interval,  implying  a  constant  radar  cross-sectional  area.  The  results  ob¬ 
tained  here  can  be  extended  to  the  fluctuating  case  using  the  Swerling  models  at  the  expense  of 
increased  complexity. 

1.2.2  Noise  Environment 

The  statistical  characterization  of  the  received  signal  is  now  examined.  The  noise  environment 
is  assumed  to  be  composed  of  Gaussian  background  noise  and/or  directional  interference  such  as 
jamming.  Therefore,  the  individual  matched  filter  outputs  are  independent  complex  Gaussian 
random  variables  (r.v.’s).  A  complex  Gaussian  r.v.  possesses  a  circularly  symmetric  pdf  with  the 
real  and  imaginary  parts  being  independent  and  each  part  having  a  variance  equal  to  (r^/2 .  Thus 
a  complex  r.v.  has  variance  cr^  given  the  preceding  definition.  The  N-element  primary  data  vector 
Zi  is  defined  to  be  complex  multivariate  Gaussian-distributed  with  its  N  x  N  covariance  matrix 
denoted  as  Mi  =  E{(2i  —  2i)(Zi  —  Zi)^}[5].  The  expected  value  of  Zj  is  0  under  H©  and  6p  under 
Hi  for  all  i.  The  term  b  is  the  complex  matched  filter  amplitude  given  that  a  target  is  present  in 
the  range  bin,  and  p  is  the  steering  vector.  As  stated  previously,  the  true  value  of  b  is  unknown. 
The  true  covariance  matrix  is  also  unknown  and  possibly  time  varying.  Hence,  the  notation  Mi 
indicates  the  time-varying  nature  of  the  noise  environment. 

The  time-varying  nature  of  the  second-order  statistic  implies  nonstationarity.  As  a  result, 
techniques  that  are  dependent  upon  wide-sense  stationarlty  assumptions  cannot  be  applied  [6].  In 
particular,  ensemble  averaging  is  not  a  valid  technique  if  the  noise  environment  is  nonstationary. 
A  nonstationary  environment  can  be  best  illustrated  by  an  example.  The  power  associated  with  a 
directional  interference  source  can  fluctuate  throughout  the  observation  Interval,  and  averages  of 
the  sample  covariance  matrix  would  not  approach  the  true  covariance  matrix  in  the  limit.  Only 
over  sampling  intervals  where  the  source  intensity  is  constant  will  an  estimate  based  upon  averaging 


^The  matched  filter  output  has  In-phase  and  quadrature  components  that  give  rise  to  the  complex 
representation. 
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be  valid.  An  extreme  case  of  a  nonstationary  environment  is  that  of  a  blinking  jammer  where  the 
jamming  power  is  randomly  switched  ON  and  OFF.  (This  scenario  is  investigated  in  Section  4.) 

In  the  sequel  no  assumptions  are  made  about  the  structure  of  the  covariance  matrix.  By 
not  Imposing  any  predefined  covariance  structure,  e.g.,  the  sum  of  dyads  plus  the  identity  matrix, 
an  unconstrained  maximization  of  the  likelihood  function  can  be  performed  with  respect  to  this 
unknown  covariance  matrix.  The  lack  of  constraints  will  be  shown  to  lead  to  the  sample  covariance 
matrix  as  the  estimate  of  the  true  covariance  matrix  at  time  instant  t. 

1.2.3  The  Decision  Rule 

Suppose  a  binary  hypothesis  test  is  to  be  performed  on  a  “data”  vector  x.  The  general  form 
of  the  pdf  of  X  is  known  but  one  or  more  of  the  true  parameter  values  of  the  distribution  are 
unknown.  To  determine  whether  the  data  belong  to  Ho  or  Hi  ,  they  are  processed  and  the  resulting 
output  is  compared  to  a  threshold.  The  GLR  test  procedure  is  a  reasonable  solution  to  a  problem 
of  this  form.  Although  not  guaranteed  to  be  optimal,  the  GLR  test  provides  a  general  formulation 
for  problems  that  contain  unknown  parameters  [2j. 

The  GLR  is  defined  as  the  quotient  of  the  joint  pdf  of  the  data  vector  x,  assuming  hypothesis 
Hi  ;  and  the  joint  pdf  of  the  data,  assuming  hypothesis  Ho.  The  terms  6i  and  ©o  are  defined  as 
the  vectors  of  unknown  parameters  under  the  respective  hypotheses.  The  numerators  and  denomi¬ 
nators  of  the  ratio  (the  likelihood  functions)  are  then  maximized  independently  with  respect  to  the 
unknown  parameters: 


GLR=  — 

0Up 

00 


p,(x;e,|Hi) 


po(x  ;  eo|Ho) 


(2) 


The  maximum  likelihood  estimates  (MLEs)  are  the  terms  that  maximize  the  respective  joint  pdf’s 
in  Elquation  (2).  For  example,  the  MLE  of  the  mean  assuming  a  Gaussian  pdf  is  simply  the  sample 
mean.  The  MLE  of  each  unknown  parameter  is  then  substituted  into  the  likelihood  ratio  in  place 
of  the  true  parameter  value.  Thus  a  decision  rule  is  obtained.  Using  x  as  the  input,  the  output 
of  the  decision  rule  is  compared  to  a  threshold.  If  the  output  exceeds  the  threshold,  Hiis  chosen, 
otherwise  the  null  hypothesis  is  chosen. 

The  statistical  performance  of  the  GLR  test  is  measured  in  terms  of  Pp  and  Pfa;  in  many 
cases  it  is  quite  difficult  to  obtain  their  exact  expressions,  and  one  may  resort  to  approximations 
or  simulations  to  obtain  these  measures  of  performance. 

1.3  Csise  Studies 

This  section  previews  the  four  cases  to  be  examined.  For  each  case,  a  brief  overview  is 
presented  and  pertinent  assumptions  or  approximations  are  stated.  The  general  form  of  the  solution 
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and  important  results  are  discussed.  This  information  can  be  found  in  Table  1  and  complete  details 
can  be  found  in  the  appropriate  sections. 


TABLE  1 


Case  Study  Summary 


Case 

No. 

of 

Pulses 

Noise 

Environment 

Comments 

IBI 

■■ 

N/A 

Pq  and  Pp^  expressed  in  terms  of  closed-form  expressions 

B 

Stationary 

Can  be  transformed  to  Case  1  via  unitary  transformation 

B 

Nonstationary 

Noncoherent  structure,  exact  Pp/^  expression,  Pq  via  simulation 

H 

■■ 

Nonstationary 

Coherent  structure.  Pq  and  Pp^  via  simulation 

1.3.1  Case  1:  The  Single  Pulse 

This  section  examines  the  GLR  test-based  decision  rule  for  one  data  sample.  This  case  has 
been  extensively  studied  by  Kelly  (7,8].  These  results  will  be  completely  reproduced  in  Section  2 
to  provide  a  framework  for  subsequent  analyses.  Exact  expressions  for  the  Pp  and  the  Pfa  are 
obtained  for  this  case.  The  Pp  is  conveniently  expressed  as  a  function  of  signal-to-noise  ratio  (SNR) 
and  the  threshold  value  corresponding  to  a  specified  Pfa  is  evaluated  using  the  incomplete  gamma 
function  [9|. 

Some  basic  terminology  is  now  defined.  A  sample  of  the  N  matched  filter  outputs  correspond¬ 
ing  to  single  return  is  denoted  as  the  N-vector  z.  The  expected  value  of  the  complex  multivariate 
Gaussian  r.v.  z  under  H]  is  6p  and  0  under  Ho.  The  complex  returned  signal  amplitude  6  is 
unknown.  In  addition,  the  covariance  matrix  under  each  hypothesis  is  unknown.  The  true  covari¬ 
ance  matrices  are  designated  as  Mi  and  Mo,  respectively.  These  unknown  parameters  must  be 
estimated  according  to  GLR  test  methodology. 

From  basic  detection  and  estimation  theory,  the  MLE  of  an  unknown  parameter  asymptot¬ 
ically  approaches  the  true  parameter  value,  i.e.,  it  is  a  consistent  estimator.  For  the  multivariate 
Gaussian  distribution,  the  MLE  of  the  covariance  matrix  is  the  sample  covariance  matrix  that  is 
a  consistent  estimator.  In  many  cases,  the  decision  rule  obtained  from  GLR  test  procedures  is 
identical  in  form  to  a  decision  rule  obtained  Wi.en  the  covariance  matrix  is  assumed  knowm  except 
that  the  sample  covariance  matrix  replaces  the  true  covariance  matrix  in  the  decision  rule.  There¬ 
fore,  the  “quality”  of  the  covariance  estimate  will  play  a  key  role  in  detector  performance  in  that  a 
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“better”  estimate  will  enable  the  GLR-based  detector  performance  to  approach  that  of  the  known 
covariance  case. 

In  an  attempt  to  improve  detector  performance  in  cases  where  the  true  covariance  matrix  is 
unknown,  Kelly  proposed  augmenting  the  primary  data  vector  with  K  additional  data  samples  that 
are  assumed  to  be  free  of  signal  returns.  These  additional  samples  are  taken  from  K  range  gates 
that  are  in  close  proximity  of  the  range  gate  currently  under  observation.  These  supplementary 
vectors  are  defined  as  secondary  vectors  and  are  denoted  as  z{k) ,  A:  =  1, . . . ,  AT.  There  are  now 
K  +  1  terms  available  for  covariance  matrix  estimation  at  a  particular  time  instant. 

Given  sufficient  range  and  spatial  resolution,  it  is  not  unreasonable  to  assume  that  multiple 
targets  would  be  sparsely  distributed.  As  a  result,  each  -element  secondary  vector  is  assumed  to 
be  zero-m  an  and  possess  the  same  covariance  matrix  as  the  primary  vector,  and  each  secondary 
vector  is  statistically  independent  of  the  primary  vector  and  the  other  K  -  1  secondaries.  These 
assumptions  are  realistic  and  significantly  reduce  the  complexity  of  the  subsequent  analysis. 

The  GLR  test  for  the  singlepulse  case  is  now  formulated.  The  respective  likelihood  functions 
are  the  joint  pdf  of  the  primary  vector  and  the  K  secondary  vectors  under  under  each  hypothesis. 
Under  Hi,  the  joint  pdf  is  the  product  of  the  primary  vector  nonzero-mean  multivariate  Gaussian 
pdf  and  K  identical  zero-mean  multivariate  Gaussian  pdf’s.  Under  Ho,  the  joint  pdf  is  simply  the 
product  of  K’  +  1  zero-mean  multivariate  Gaussian  pdf’s. 

The  GLR  test  procedure  requires  that  the  unknown  parameters  be  replaced  by  their  MLEJs. 
First,  each  likelihood  function  is  maximized  with  respect  to  the  unknown  covariance  matrix  while 
holding  the  returned  signal  amplitude  fixed.  The  resulting  MLEs  are  the  sample  covariance  ma¬ 
trices.  Note  that  the  GLR  test  procedure  leads  directly  io  the  sample  covariance  matrix  form, 
which  is  an  unbiased  estimator  and  asymptotically  approaches  the  true  covariance.  Therefore,  the 
inclusion  of  the  K  additional  secondary  vectors  yields  a  more  accurate  estimate  of  the  covariance. 

Estimates  of  the  covariance  matrices  are  substituted  into  the  likelihood  functions  in  place 
of  the  true  covariance  matrices.  Given  the  multivariate  Gaussian  distribution  assumption,  the 
likelihood  ratio  at  this  point  is  the  quotient  of  the  determinant  of  the  MLE  of  the  covariance  under 
Ho  and  the  determinant  of  the  MLE  of  the  covariance  matrix  under  Hj.  The  MLE  of  the  covariance 
matrix  under  Hi  is  still  a  function  of  the  unknown  returned  signal  amplitude.  This  problem  is  now 
equivalent  to  the  Wilks’  Lambda  test  commonly  found  in  statistical  literature  [10].  For  this  case, 
the  MLE  of  b  is  equivalent  in  form  to  the  minimum  variance  unbiased  estimator  of  the  returned 
amplitude.  Substituting  the  MLE  of  b  into  the  test  yields  the  test  form  that  will  be  the  basis  of 
the  single-pulse  analysis. 

Under  Hi,  the  final  form  of  the  single-pulse  test  is  statistically  equivalent  to  the  inverse  of 
a  noncentral  0  r.v.  conditioned  on  the  /^-distributed  loss  factor  p.  This  loss  factor  arises  from 
the  need  to  estimate  the  covariance  matrix.  Reed,  Mallet,  and  Brennan  have  extensively  studied 
this  r.v.  and  have  shown  that  estimating  the  covarlr.nce  matrix  effectively  reduces  the  available 
SNR  to  the  detector,  resulting  in  performance  degradation  (11).  Section  2  shows  that  as  the 
number  of  secondary  vectors  used  for  covariance  estimation  increases,  estimation  loss  decreases 
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and  detector  performance  approaches  that  of  the  known  covariance  (matched  filter)  case.  Under 
the  null  hypothesis,  the  test  is  statistically  equivalent  to  a  central  0  r.v.  Also,  exact  expressions 
for  the  PpA  and  Pd  are  derived  in  Section  2. 

1.3.2  Case  2:  Multiple  Pulses  in  a  Stationary  Environment 

The  test  for  this  case  uses  Np  observation  Intervals  to  make  a  decision.  By  making  the 
following  assumptions,  this  problem  can  be  transformed  into  a  structure  identical  to  that  described 
in  case  1.  First,  the  complex  returned  amplitude  b  is  assumed  constant  for  the  duration  of  the 
test.  Second,  for  this  case  a  wide-sense  stationary  noise  environment  is  a.ssumed  over  the  entire 
observation  interval;  therefore,  the  second-order  noise  statistic,  the  true  covariance,  does  not  vary 
as  a  function  of  time.  The  data  are  assumed  to  be  statistically  independent  of  the  noise  and 
interference. 

Suppose  an  N  y.  Np  matrix  is  formed,  consisting  of  Np  statistically  independent  primary 
vectors  and  collect  K  secondary  data  vectors  during  the  observation  intervals.  It  can  be  shown 
that  postmultiplication  of  this  primary  data  matrix  by  an  appropriate  data-dependent  unitary 
transformation  results  in  the  rotation  of  the  mean  value  component  of  each  Np  primary  vector  into 
the  primary  data  matrix  first  column  position  (12).  The  remaining  Np  -  1  primary  vectors  are  now 
signal-free  or  equivalently,  they  are  zero-mean  vectors.  Thus  the  Nppulse  test  has  been  transformed 
into  the  single-pulse  test  form  but  with  K  +  Np  —  I  zero-mean  secondary  vectors.  The  same  Pd 
and  PpA  formulas  as  derived  for  case  1  can  be  applied. 

For  large  Np,  the  additional  Np  -  1  vectors  used  for  covariance  estimation  will  significantly 
reduce  estimation  loss.  Of  course,  secondaries  obtained  from  each  observation  Interval  (NpxK)  can 
be  used  for  the  covariance  estimate,  but  the  additional  memory  and  computational  requirements 
could  be  expensive.  The  statlonarlty  assumption  allows  the  primary  data  vectors  themselves  to  be 
used  to  Improve  the  covariance  estimate. 

Because  the  case  2  test  can  be  transformed  into  a  previously  solved  problem,  the  analysis  is 
Included  In  Section  2.  In  either  case,  improved  covariance  estimation  results  in  improved  detector 
performance. 

1.3.3  Multiple  Pulses  in  a  Nonstationary  Environment 

When  considering  the  multiple-pulse  test,  a  nonstationary  noise  environment  rules  out  the 
solution  technique  used  for  case  2.  A  unitary  transformation  cannot  be  applied  to  the  primary  data 
matrix  because  each  vector  is  no  longer  statistically  Identical;  therefore,  the  discussion  continues 
with  the  general  formulation  of  the  GLR  for  multiple-pulse  conditions  where  there  ate  numerous 
unknown  parameters. 

The  sampled  data  are  assumed  statistically  independent.  Due  to  this  assumption,  the  like¬ 
lihood  functions  used  to  form  the  GLR  are  simply  the  product  of  the  Np  likelihood  functions 
associated  with  each  data  sample.  If  the  noise  is  Gaussian,  each  of  the  Np  likelihood  functions 


7 


takes  on  the  form  of  the  multivariate  Gaussian  distribution.  GLR  test  procedures  indicate  that 
the  likelihood  function  for  the  complete  data  set  under  each  hypothesis  must  then  be  maximized 
with  respect  to  each  unknown  parameter.  In  the  case  where  each  unknown  parameter  varies  as  a 
function  of  time,  maximization  of  the  overall  likelihood  function  is  equivalent  to  maximizing  each 
of  the  Nj,  likelihood  functions  independently.  If  one  or  more  of  the  unknown  parameters  remains 
coastant,  however,  the  maximization  process  can  become  quite  difficult  or  nearly  impossible. 

For  the  problem  under  investigation,  the  true  covariance  matrix  is  assumed  to  be  unknown 
and  time  varying.  As  previously  discussed,  the  signal  amplitude  is  constant  over  the  complete  ob¬ 
servation  interval.  The  MLE  of  each  of  the  Np  covariance  matrices  can  be  evaluated  by  maximizing 
the  likelihood  function  under  each  hypothesis  at  each  time  instant.  The  MLE  of  each  covariance 
matrix  is  the  sample  covariance  matrix  obtained  iLsing  the  primary  vector  and  the  K  secondary 
vectors;  however,  because  the  signal  amplitude  is  constant,  the  maximization  of  the  GLR  with 
respect  to  this  unknown  parameter  requires  the  solution  of  a  2A/p-degree  polynomial.  For  even  a 
relatively  small  Np,  the  problem  is  unyielding. 

Two  possible  techniques  were  developed  to  deal  with  this  problem.  First,  allow  the  model  for 
the  returned  amplitudes  be  time  varying.  This  modeling  assumption  enables  the  likelihood  ratio 
to  be  expressed  as  the  product  of  Np  independent  ratios,  which  results  in  the  same  form  as  the 
ratio  found  in  Equation  (2).  This  reformulation  of  the  problem  is  the  basis  of  the  case  3  analysis. 
Second,  do  not  allow  the  returned  amplitude  model  to  be  time  varying  (the  actual  condition  for  a 
constant  radar  cross-sectional  area).  After  maximizing  the  test  with  respect  to  the  Np  unknown 
covariances,  the  test  can  be  expressed  in  terms  of  a  series  expansion  using  an  algebraic  identity. 
The  series  expansion  approach  is  the  basis  of  the  case  4  analysis.  These  two  cases  are  now  discussed 
in  greater  detail. 

Case  S  :  Noncoherent  Detection  Using  Multiple  Pulses.  By  allowing  the  model  for  the  signal 
mplltude  to  be  time  varying,  the  decision  rule  can  be  expressed  as  the  product  of  Np  of  the 
livelihood  ratios  described  in  case  1  [13].  This  assumption  enables  independent  maximization  of 
ei.ch  Np  likelihood  function  with  respect  to  the  unknown  parameters.  The  GLR  for  this  case  is 
sti  listically  equivalent  to  the  product  of  Np  Independent  /S-distributed  r.v.’s.  Section  3  will  show 
thar  the  test  is  equivalent  to  a  noncoherent  integration  process.  Performance  of  this  noncoherent 
detector  is  expected  to  approach  the  matched  filter  performance  level  as  K  -*  oc. 

To  evaluate  noncoherent  detector  performance,  the  threshold  value  corresponding  to  a  speci¬ 
fied  PpA  must  be  determined.  Fortunately,  using  the  transformation  of  variable  technique  [6],  the 
product  of  Independent  central  (3  r.v.’s  with  a  special  set  of  degrees-of-freedom  parameters  can  be 
transformed  into  a  central  degrees  of  freedom.  The  threshold  value  is  obtained 

using  the  incomplete  gamma  function  (14). 

The  Pd  for  this  test  presents  considerable  difficulties.  For  Np  >  2,  it  is  virtually  Impos¬ 
sible  to  obtain  an  exact  expression  for  the  distribution  of  the  product  of  noncentral  0  r.v.’s.  A 
Chemoff  bound  for  this  distribution  was  obtained  and  provided  an  analytical  means  for  evaluat¬ 
ing  performance.  This  approach  is  quite  different  from  the  approximate  x*  techniques  preferred  by 
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statisticians  when  analy2ing  Wilks’  Lambda  tests  [15] .  A  Monte  Carlo  simulation  was  used  to  obtain 
experimental  results.  As  expected,  the  analysis  revealed  that  noncoherent  detector  performance 
Improved  as  a  function  of  Np  and  K. 

Case  4-‘  Coherent  Detection  Using  Multiple  Pulses.  Unlike  case  3,  the  model  for  the  returned 
signal  amplitude  is  held  constant  over  the  observation  interval.  The  GLR  test,  as  per  case  3,  can  be 
maximized  with  respect  to  the  Np  unknown  covariance  matrices;  however,  the  maximization  of  the 
GLR  with  respect  to  b  is  almost  impossible.  This  problem  can  be  solved  using  an  approximation. 
The  GLR  is  expressed  in  terms  of  a  truncated  power  series  expansion.  The  approximate  GLR  can 
then  be  maximized  with  respect  to  6.  After  some  manipulation,  the  final  approximate  GLR  can 
be  expressed  as  the  product  of  a  term  that  is  a  function  of  the  amplitude  estimate  and  the  GLR 
obtained  for  case  3.  This  result  is  defined  as  a  coherent  decision  rule  because  the  estimate  of  b 
incorporates  data  from  the  complete  Np  pulse  observation  interval.  The  statistical  analysis  of  this 
decision  rule,  however,  remains  daunting  and  analytical  attempts  were  abandoned  in  favor  of  Monte 
Carlo  simulations. 

Several  conditions  were  examined  via  simulation.  As  a  performance  benchmark,  the  coherent 
detector  was  tested  in  the  less  severe  stationary  noise  environment.  Given  a  reasonable  number 
of  secondary  vectors,  detector  performance  exceeded  that  of  the  noncoherent,  known  covariance 
detector.  Under  nonstationary  conditions,  the  coherent  detector  also  exceeded  the  performance  of 
the  noncoherent  detector.  In  fact,  coherent  detector  performance  decreased  only  when  the  number 
of  covariance  estimates  during  an  Np  pulse  observation  interval  was  reduced.  Among  the  noise 
environments  used  for  test  purposes  were  the  blinking  jammer  and  multiple  jammer  scenarios. 

1.4  Report  Outline 

The  remainder  of  this  report  is  organized  as  follows;  Section  2  presents  the  single-pulse 
analysis  based  upon  the  GLR  test.  Because  the  multiple-pulse  stationary  case  can  be  reduced 
to  the  form  of  case  1,  its  analysis  is  Incorporated  here.  Section  3  explc  s  the  performance  of  the 
noncoherent  detector,  including  a  complete  derivation  of  the  Chernoff  bound.  The  coherent  detector 
that  is  based  upon  a  series  expansion  is  examined  in  Section  4.  The  final  section  summarizes  the 
results  and  conclusions  for  all  cases  (see  Table  1). 
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2.  SINGLB-PULSE  ADAPTIVE  DETECTOR 


In  the  course  of  deriving  the  GLR  test,  basic  terminology  is  presented.  The  resulting  test  is 
shown  to  be  statistically  equivalent  to  1  plus  a  conditional  F-distributed  r.v.  Using  a  relationship 
between  F  and  ^  r.v.’s,  the  inverse  of  the  test  statistic  is  shown  to  possess  a  conditional  /?  distribu¬ 
tion.  U.'-ing  the  cumulative  distribution  function  of  the  0  distribution,  exact  expressions  for  the  Pd 
and  PpA  are  obtained,  and  the  effect  of  covariance  estimation  on  detector  performance  is  discussed 

17]. 


2.1  Single-Pulse  GLR  Test 

The  analysis  begins  by  defining  the  complex  multivariate  Gaussian  pdf  that  characterizes  the 
primary  and  secondary  data  vectors.  As  stated  in  Section  1,  the  real  and  imaginary  parts  of  each 
element  of  the  data  vectors  are  independent.  Under  the  signal  plus  noise  hypothesis,  the  mean  of 
the  primary  data  vector  is  bp;  and  under  the  noise-only  hypothesis  the  mean  vector  is  0.  The  pdf 
of  the  primary  vector  under  Hi  is 


where 


Mi  =  {E(z-hp)(2-bp)«}  , 

b  Is  the  complex  returned  signal  amplitude,  and  p  is  the  steering  vector.  Without  loss  of  generality, 
the  steering  vector  is  defined  to  have  unit  length. 

The  primary  vector  pdf  under  Ho  Is 


Po(2lHo)  = 


(4) 


where 


Mo  =  e{zz^} 

The  secondary  vectors,  zik),k  =  are  also  zero-mean  and  possess  the  same  pdf  as  the 

primary  vector  under  Ho;  Mi  and  Mo  are  assumed  unknown. 
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Prom  previous  studies  where  covariance  matrices  are  assumed  unknown,  detector  performance 
was  shown  to  be  directly  related  to  the  quality  ol  the  covariance  estimate  (7,11).  Because  the  sample 
covariance  matrix  is  the  estimator  form  obtained  from  the  GLR  test  procedure,  this  estimate  will 
asymptotically  approach  the  true  value  as  the  number  of  sample  vectors  used  to  form  the  estimate 
increases.  Therefore,  the  covariance  matrix  estimator  should  use  the  maximum  number  of  available 
samples. 

The  joint  pdf  of  the  primary  and  K  secondary  vectors  forms  the  basis  for  the  GLR  procedure. 
Because  nearby  range  gates  are  assumed  to  share  the  same  covariance  matrbc,  K  +  1  data  vectors 
possessing  the  same  covariance  are  now  available  to  form  a  decision  rule.  The  primary  and  secondary 
vectors  are  statistically  independent;  hence,  the  joint  pdf  of  the  primary  vector  and  the  K  secondary 
vectors  is  simply  the  product  of  the  Individual  pdf’s.  The  joint  pdf  is  denoted  as 

Pi{z,2(l),...,z(A:)|Ha,i  =  0,l  .  (5) 


Forming  the  ratio  of  the  respective  joint  pdf’s,  the  single-pulse  likelihood  ratio  is 


Pi{»,2(1),...,z(A-)|Hi} 

Po{z,2(1),...,z(A:)|Ho} 


The  numerator  and  denominator  pdf’s  are  referred  to  as  the  “likelihood  functions.” 

Equation  (6)  is  a  function  of  the  unknown  parameters  6,  Mi,  and  Mo-  Because  of  the  presence 
of  unknown  parameters,  the  standard  GLR  test  procedure  Is  used  to  arrive  at  a  decision  rule.  The 
GLR  expressed  In  terms  of  the  multivariate  Gaussian  pdf  Is  defined  as: 


A(M,,Mo,fe)  = 


•up 

Ml, 6 


flUp 

Mo 


^-{z"Mo*z+Q} 


(7) 


where 


Ci  = 


1 


IIMill 


,i  =  0,l 


and 
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Q 


K 


Y^zf^ik)M^^z{k) 


For  mathematical  simplicity,  apply  the  inner  product  identity  x^Y"*x  =  Tr{Y~*xx^}  to  the 
{K  +  l)th  root  of  Equation  (7).  Using  this  identity  yields 

aup  Cl  e  I  i 

-TV 

8up  Co  e 
Ml 


where 


■  S  =  '^z{k)z^ik) 

isl 


A  is  now  redefined  as  the  {K  +  l)th  root  of  the  ratio. 

The  GLR  test  procedure  requires  that  the  MLEs  of  the  unknown  parameter  under  each 
hypothesis  be  evaluated  and  then  substituted  Into  the  test  in  place  of  the  true  parameter  values. 
The  preceding  expression  is  solved  to  obtain  the  MLEs  of  the  unknown  covariance  matrix.  The 
MLEs  are  denoted  as  Mi  and  Moi  respectively.  From  Equation  (8),  the  MLEs  of  the  covariance 
matrix  under  each  hypothesis  are 

M,(6)  =  {(z  -  6p)(z  -  6p)«  +  S  }  (9) 


and 


Mo  =  -j^{xi''+s}  (10) 

As  stated  in  Section  1,  the  sample  covariance  form  found  in  Equations  (9)  and  (10)  arises  from  the 
problem  formulation,  not  from  a  priori  assumption. 

The  MLE  Mi  is  still  a  function  of  the  unknown  returned  amplitude  that  must  be  subsequently 
estimated.  After  substituting  the  MLEs  of  the  covariance  matrices  into  their  respective  likelihood 
functions,  the  GLR  reduces  to  the  following  ratio  of  determinants: 


13 


(11) 


A(6)  =  Jfe!L 

min  II  Ml  II 

b 

The  process  of  minimizing  the  denominator  is  equivalent  to  maximizing  the  likelihood  ratio.  Using 
the  following  determinant  identities, 

(/C  +  l)^||M,||  =  ||S||(l  +  (z-fcp)"S-»(2-i>p)l  (12) 


and 


(/f  +  l)'^l|Moll  =  ||S|![l  +  2"S-*zl  ,  (13) 

and  completing  the  square  with  respect  to  b,  results  in  an  expression  where  the  unknown  parameter 
b  is  isolated.  The  GLR  can  be  written  as 


A(6)  = 


1  +z^S"*z 


mln^  J+Z^S“*Z+p^S*  *p 

b 


p«S-‘z 

1 

p"S-‘z 

Cl 

p"S-*p 

The  denominator  in  Equation  (14)  is  minimized  when  b  =  b,  where  6  is 


6  = 


p«S-‘z 

P^S-‘p 


(14) 


(15) 


The  MLE  b  is  identical  in  form  to  the  minimum  variance  unbiased  estimator  of  b  given  that  l/K 
S  replaces  M.  The  estimate  is  essentially  produced  by  whitening  and  cross-correlating  the  data 
vector  with  the  steering  vector.  The  substitution  of  b  into  the  GLR  yields 


A- 


l+z^S-'z 


l-l-z«S-‘z- 


p^S'‘z 


> 

■f  <  a 


(16) 


where  a  is  a  threshold  parameter. 

Before  deriving  the  statistics  of  the  GLR  test,  an  attempt  was  made  to  develop  an  intuitive 
feel  for  the  behavior  of  the  test.  As  found  in  Equation  (16),  the  test  is  at  a  maximum  when 
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the  denominator  term  takes  on  its  minimum  value.  The  conditions  under  the  signal-plus-noise 
hypothesis  are  examined.  By  setting  the  primary  data  vector  z  equal  its  expected  value  under  Hj, 
ftp,  and  performing  some  algebraic  manipulation,  the  denominator  reduces  to  the  scalar  value  1, 
which  is  the  minimum  value  of  the  denominator  for  any  data  vector  z.  If  the  target  is  present 
and  is  collinear  in  space  with  the  steering  vector,  the  best  chance  of  making  the  correct  decision  is 
at  hand.  The  numerator  term  under  this  assumption  takes  on  the  form  1  4-  |6|^p^S“^p.  As  the 
sample  covariance  approaches  the  true  covariance,  |6pp^S“‘p  becomes  in  the  limit  |6|*p^M~*p, 
which  is  the  nonadaptive  SNR  of  the  array  except  for  the  factor  XjK.  Therefore,  as  the  covariance 
estimate  improves  the  SNR  required  to  achieve  a  certain  Pp  decreases. 

Alternatively,  the  presence  of  the  sample  covariance  in  the  decision  rule  results  in  the  actual 
SNR  seen  by  the  detector  to  be  less  than  the  optimal  nonadaptive  SNR.  If  the  estimation  loss  factor 
is  defined  to  be  the  random  variable  p ,  the  performance  of  the  test  can  be  evaluated  in  terms  of 
the  product  of  the  nonadaptive  SNR  and  the  complex  /3-distributed  loss  factor  p.  This  loss  factor 
was  derived  in  Reed,  Mallett,  and  Brennan  [111,  and  its  effect  on  the  performance  of  the  decision 
rule  was  first  studied  in  Kelly  (7).  The  loss  factor  plays  a  critical  role  in  the  performance  of  the 
adaptive  detector. 

2.2  Statistical  Representation  of  the  Single^Pulse  GLR  Test 

The  analysis  of  the  performance  of  the  single-pulse  GLR  test  is  based  upon  several  standard 
multivariate  statistical  techniques  [16-18].  The  primary  and  secondary  data  vectors  are  whitened 
so  that  each  of  the  components  of  the  N-element  data  vectors  are  statistically  Independent.  The 
space  spanned  by  the  data  vectors  is  then  divided  into  two  orthogonal  subspaces  of  dimension  1  and 
N-l,  respectively.  A  unitary  transformation  is  then  applied  to  the  data  vectors  to  rotate  the  mean 
value  of  the  data  completely  into  the  first  element  of  the  data  vectors.  By  fixing  the  components  of 
the  N  —  1  element  subspace,  the  conditional  distribution  of  the  single-element  subspace  is  obtained. 
After  removing  the  conditioning  effect  through  averaging,  the  unconditional  distribution  of  the  test 
provides  the  basis  for  the  Pd  and  Pfa  expressions. 

For  clarity,  the  GLR  test  is  restated: 


A  = 


1  +  z^S"*z 


l+z«S-‘z- 


p'^S-* 


> 

y  <  Qt 


(17) 


and  the  analysis  proceeds  by  whitening  the  data  vectors.  An  appropriate  choice  for  a  whitening 
transformation  is  the  inverse  of  the  square  root  of  the  covariance  matrix.  Because  the  covariance 
matrix  (and  its  Inverse)  is  positive  definite,  a  positive  definite  square  root  matrix  can  be 

introduced.  By  defining  the  whitened  primary  vector  vector  as  z  =  z,  then  Ez  = 

E I z)  ~  choice  of  as  a  transformation  matrix  results  in 
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a  statistically  lndei>endent  data  set.  Likewise,  a  similar  definition  of  the  whitened  secondary  vector 
z(Jfc)  yields  the  same  result.  The  steering  vector  is  also  whitened  to  obtain  p  = 

A  unitary  transformation  Ur  is  defined  such  that  the  whitened  steering  vector  p  = 
is  transformed  into  a  basis  vector  with  1  as  the  first  coordinate  tuid  0  as  the  remaining  N  -  I 
elements.  This  rotation  can  be  accomplished  by  selecting  the  first  row  of  U/i  to  be  a  normalized 
unit  vector  and  the  other  rows  of  U/j  to  be  orthogonal  to  the  first  row  and  to  each  other,  i.e., 

=  6ij,  where  6ij  is  the  Kronecker  delta  function.  The  whitened,  rotated  steering  vector  is  defined 
as 


e  =  Ufip  =  (l  (18) 

The  test  statistic  is  unchanged  by  a  unitary'  transformation  because  it  does  not  change  a 
vector  norm.  The  premultiplication  of  z  and  i(k)  by  Ur  affects  only  the  mean  of  z.  To  simplify 
notation,  the  whitened,  rotated  primary  and  secondary  data  vectors  are  redefined  in  terms  of  the 
original  notation  as  the  vectors  z  and  2(fc),  respectively.  TheGLR  in  terms  of  the  whitened,  rotated 
vectors  is 


A  =  - 


1  +  z"  S-‘z 


1+z"  S-‘z~ 


(19) 


where  now 


S  =  Ur  S  Ufi^ 


and 


E{z}  =hUR  M-‘/2p  . 

The  data  vectors  are  partitioned  into  an  A  component  that  corresponds  to  the  first  element 
and  a  B  component  that  corresponds  to  the  remaining  N  —  1  elements.  The  partitioned  primary 
and  secondary  vectors  are 


*  =  =  Ur  Z 


(20) 


where  za  Is  ®  scalar,  zr  is  a  1  x  (N  -  1)  vector,  and 
z(A:)  =  (za(^)  =  Ur  z(/c) 
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The  choice  of  the  transformed  steering  vector  p  results  in  this  particular  partitioning.  The 
whitened,  rotated  S  matrix,  which  is  K  times  the  sample  covariance  matrix  based  on  the  secondaries 
alone,  is  partitioned  as 


where 


« _  ^AA  Sab 

O  —  ) 

Sfi/i  Sbb 


S/u  =  j3  XA{k)zA^ik)  , 


Sab  =  XI  , 

i=l 


Sba  =  ^B{k)zA^{k)  , 


Sbb  =  ^  ZB{k)zB'^{k) 


The  inverse  of  the  sample  covariance  matrix  is  defined  as  P  and  is  partitioned  in  a  manner  analogous 
to  that  used  for  the  matrix  S,  l.e., 


Paa  Pab 


Pba  Pbb 


The  three  inner  product  terms  of  the  GLR  test  are  now  expressed  in  terms  of  the  preceding 
definitions.  First, 

z"S-‘2  =  z^Pz=  (23) 


(*A  +  PaaP/ib^o)  Paa  (*a+PaaPabzb) 

+  *B  (PbB  -  PbaPaaPab)  zb 


Second, 
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e«S-‘e  =  e"Pe  =  PAyi  , 


(24) 


And  third,  the  final  inner  product 

e^S-*z  =  e^Pz  =  P^,4  (2,4  +  P^Pabzb)  (25) 

Using  the  Schur  relations  for  partitioned  matrices,  the  inner  product  z^S'^z  is  redefined: 

z«S-»z  =  zgSi^ZB  +  ^  ,  (26) 

where 

V  =  ^A-  SxfiSgjjZfi  , 

and 

T  =  SxM  —  SabSqqSba  —  Paa 

It  is  interesting  to  note  that  terms  y  and  b  can  be  shown  to  be  identical  using  Ekiuations  (26) 
and  (27). 

To  digress  slightly,  the  discussion  turns  to  the  significance  of  the  r.v.’s  y  and  T.  Suppose 
that  a  Gaussian-distributed  N-vector  g  has  a  mean  /x  and  covariance  matrix  E.  The  vector  g  is 
partitioned  into  q  and  N  —  q  dimensional  subvectors  and  the  covariance  niatrix  E  is  partitioned  in 
the  same  way  as  S.  A  linear  transformation  to  the  new  components  hi  and  h2,  which  results  in  hi 
being  uncorrelated  with  h2  is 

hi  =  gi  -  Ei2E^‘g2  (27) 

and 

h2  =  g2  ,  (28) 

where  E  is  partitioned  as 

E= 

E2I  E22 
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The  product  S12  Sjj'  is  referred  to  as  the  “regression  coefficient”  in  the  statistical  literature 
[16j.  From  a  filtering  or  estimation  theory  viewpoint,  hi  is  the  prediction  error  given  h2  . 

The  conditional  distribution  p(hi  |h2)  is  Gaussian-distributed  with  mean  n\  - 
and  covariance  En  -  Ei2Ej2*E2i.  Observing  that  the  form  of  y  is  identical  to  that  of  hi,  the  r,v. 
y  is  simply  the  prediction  error  of  za  given  zb-  The  conditional  covariance  is  not  a  function  of 
h2;  therefore,  En  -  E12  E22  S21  is  the  unconditional  covariance  of  hi-  Because  T  is  Identical 
in  structure  to  the  unconditional  covariance  of  y,  T  is  clearly  independent  of  y.  The  preceding 
discussion  can  be  used  to  gain  valuable  insight  into  the  final  form  of  the  decision  rule. 

The  remaining  inner  products  expressed  in  terms  of  the  r.v.’s  y  and  T  are 

e«S-*e  =  e«Pe=  i  (29) 

and 

e"S-4=eWpz  =  i  (30) 

The  GLR  expressed  in  terms  of  the  Schur  components  is 
1  +  2^S~*Z 

A  =  --  -g-n —  •  (31) 

1  +  XgSggZB 


Equation  (31)  is  now  used  along  with  the  definitions  of  y  and  T  to  obtain  the  final  form  of 
the  GLR.  The  r.v,  v  is  defined  as 


V  = 


[l+zgS^^ZB 
where  the  loss  factor  p  is  defined  as 
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(32) 


[1  +  Zb^^^Zb] 

The  GLR  test  for  the  single-pulse  case  expressed  in  terms  of  t;  and  T  is 


(33) 


A  =  l4-i^^a  . 


(34) 
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Equation  (32)  is  the  form  of  the  decision  rule  used  to  obtain  test  performance  criteria,  i.e., 
Pd  and  Pfa- 

The  statistical  characterization  of  the  previous  test  is  now  examined.  In  general  terms,  given 
a  vector  x  of  length  N  distributed  as  and  a  sample  covariance  matrix  S  composed  of  K 

zero-mean  terms,  the  quadratic  inner  product  x^S“*x  is  distributed  as  noncentral  F  with  N  and 
K  -  N  +  1  complex  degrees  of  freedom  (dof)  and  noncentrality  parameter  [16].*  This 

definition  of  the  complex  Gaussian  r.v.  accounts  for  the  correspondence  of  one  complex  to  two 
real  dof;  therefore,  the  pdf  for  a  real  r.v.  with  an  even  number  of  dof  is  equivalent  to  the  pdf 
of  complex  r.v.  with  one-half  the  number  of  dof.  Fbr  distributions  based  upon  the  multivariate 
Gaussian  distribution,  the  noncentrality  parameter  is  the  squared  norm  of  the  mean  vector.  If  x  is 
zero-mean,  the  resulting  distribution  is  central  F  with  the  same  dof.  The  noncentrality  parameter 
is,  of  course,  0  for  central  distributions.  Therefore,  Zg  SggZB,  the  denominator  of  p,  is  distributed 
as  central  F  with  N  —  I  and  K  —  N  +  2  complex  dof.  An  important '  lationship  between  F  and  /? 
r.v.’s  is 


xp(n;  m) 


1 

1  -i-XF(m,n) 


(35) 


where  m  and  n  are  the  dof  parameters;  therefore,  the  loss  factor  p  is  distributed  as  central  0  with 
K  -N  +  2&x\d  N  ~l  dof. 

Because  v  is  a  conditional  Gaussian  r.v.,  the  likelihood  ratio  term  |v|^  is  ^^-distributed  con¬ 
ditioned  on  the  /?-distributed  loss  factor  p-  Because  v  is  composed  of  only  a  single  term,  it  has 
just  one  complex  dof.  In  terms  of  the  original  coordinates,  the  mean  vector  of  z  is  The 

squared  norm  of  this  vector  is  (6M“*''^p)^(6M“‘/V)  =  which  is  simply  the  non- 

adaptlve  SNR.  Under  conditioning,  the  B  components  are  held  constant  and  the  conditional  mean 
of  V  is  the  product  bM~U2ppi/2_  resulting  noncentrality  parameter  under  conditioning  is  the 
product  of  the  nonadaptive  SNR  and  the  loss  factor  p.  The  loss  factor  enters  the  test  statistic  only 
through  the  conditional  mean. 

The  r.v.  T,  which  equals  Saa  -  SabSb^Sba.  Is  distributed  as  central  with  K  ~  N  +  1 
dof  (5|.  The  distribution  of  T  can  be  obtained  using  the  prop>erties  of  Wishart  matrices  and  certain 
orthogonal  transformations.  T  is  not  a  function  of  the  primary  vector  subcomponents,  hence  it  is 
statistically  independent  of  y  and  v. 

The  ratio  Irj^/T  is  composed  of  two  Independent  r.v.’s.  Because  the  ratio  of  independent 
r.v.’s  is  defined  to  be  F-distrlbuted,  the  term  |v|^/r  is  conditionally  F-distributed  with  1  and 
K  —  N  +  1  dof.  Using  the  relationship  expressed  in  Equation  (35),  A  is  statistically  equivalent 


^Derivations  are  for  real  r.v.’s  but  can  be  extended  to  the  complex  r.v.  case. 
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to  the  inverse  of  a  conditional  y9-distributed  random  variable.  The  test  statistic  is  now  examined 
under  hypotheses  Hi  and  Ho- 

2.3  Single-Pulse  Pp;^  and  Pq 

To  obtain  expressions  for  the  Pp  and  Pfai  the  conditional  and  unconditional  distributions  of 
the  test  statistic  are  examined  first.  As  stated  in  Section  2.2,  the  inverse  test  statistic  is  complex 
/3*di8tributed.  The  conditional  distribution  of  the  inverse  test  statistic  A“*  is 

P(x1p)  =  f/3(a:;L,  l;c)  ,  (36) 


where 


L=K-N+1  , 

c  =  ap=  p  under  Hi 


and 


c  =  0  under  Ho 

The  presence  of  the  conditioning  r.v.  requires  that  the  distribution  of  the  test  statistic,  which 
is  conditionally  /?-dlsirlbuted,  is  averaged  over  p  to  obtain  the  unconditional  distribution  of  the 
test.  Under  Ho  the  noncentrality  parameter  is  0;  therefore,  the  conditioning  r.v.  has  no  effect  on 
the  unconditional  distribution  because  it  appears  only  through  the  noncentrality  parameter.  As  a 
result,  the  inverse  test  distribution  is  unconditionally  central  /^-distributed  under  Ho. 

Because  its  form  is  much  simpler,  the  distribution  of  the  test  under  Ho  is  examined  first.  The 
analysis  results  in  an  exact  expression  for  the  Pra-  The  general  form  of  the  central  /?  distribution 
is  [5] 


fp(a:;  m,n) 


r(n-t-m)  ^1 
r(n)r(m)"^ 


(i-i)"'' 


(37) 


where  n  and  m  are  the  complex  dof  parameters. 

In  the  current  case  where  the  first  dof  parameter  is  L  and  the  second  is  1,  the  central  /? 
distribution  reduces  to  the  simple  form 


fp(x;L,  1)  = 


(38) 
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The  Pfa  is  defined  as: 


Pfa=  /~f(A|Ho)dA  .  (39) 

Ja 

Using  the  relationship 

Pr(A  >a)  =  Pr(A-*  <  ,  (40) 

and  the  fact  that  A~*  is  unconditionally  distributed  as  {0{x;  L,  1),  Pfa  is  the  cumulative  distribution 
function  of  a  central  /?  r.v.,  i.e., 


Pfa  =  r  Lx^-^dx  =  .  (41) 

Jo 

Because  Pfa  is  independent  of  the  true  covariance  matrix  and  is  only  a  function  of  the  dimensional 
parameters  of  the  problem,  the  test  statistic  possesses  the  desirable  constant  false  alarm  rate 
(CFAR)  characteristic.  The  threshold  value  a  is  evaluated  solving  Equation  (41)  for  a  sF)ecific  Pfa 
value. 

Under  the  signal-plus-noise  hypothesis,  the  presence  of  the  conditioning  r.v.  results  in  a  more 
complicated  analysis.  The  unconditional  distribution  of  the  test  statistic,  which  is  obtained  by 
averaging  out  over  the  distribution  of  the  loss  factor  (contained  in  c),  is 

f(i|H,)=  f\0{x-,L,l,c){0{p-,L  +  l,N-l)dp  .  (42) 

Jo 

Per  the  Pfa  calculation,  the  Pd  can  be  expressed  in  terms  of  the  unconditional  test  distribu¬ 
tion  as 


-1 


Pd  =  Pr(A  >a)=f  f(x|Hi)dx 
Jo 


Using  previously  derived  results  (8),  Pd  Is 


where  a  is  the  nonadaptive  SNR, 


(43) 


(44) 
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and 


Hfc(u;)  =  /'  Gfc(ap)  (p;  I  +  1,  AT  _  i)  dp 
Jo 


k-\ 


Gt(l)  =  e- 

m\ 


maO 


A  complete  discussion  of  the  computational  complexities  involved  in  the  evaluation  of  the  proba¬ 
bility  of  detection  can  be  found  elsewhere  (19,20). 

As  discussed  in  greater  detail  in  the  remainder  of  the  report,  the  formulation  of  the  GLR 
test  is  equivalent  to  the  Wilks’  A  test,  which  commonly  occurs  in  multivariate  hypothesis  testing 
schemes  (10).  Results  obtained  for  the  signal-plus-noise  hypothesis  differ  from  the  standard  results 
because  of  the  presence  of  the  conditioning  r.v.  Next  is  an  examination  of  the  performance  of  the 
decision  rule  as  a  function  of  the  nonadaptive  SNR  and  the  excess  number  of  secondary  vectors  L. 

2.4  Performance  of  the  Single-Ptilse  Detector 

The  statistical  performance  of  the  single-pulse  decision  rule  is  quantified  through  the  use  of 
receiver  operating  characteristic  (ROC)  plots.  ROC  is  defined  here  as  the  plot  of  the  Pd  versus 
SNR  for  a  fixed  Pfa,  which  is  set  to  a  nominal  value  of  10“^.  The  parameter  L,  which  is  related 
to  the  number  of  secondary  vectors,  is  a  controlling  factor  in  the  performance  of  the  decision  rule. 
As  L  tends  to  oo,  the  loss  factor  tends  to  1,  hence,  the  covariance  estimate  approaches  the  true 
covariance.  When  there  is  no  estimation  loss,  the  Marcum-Q  function  can  be  used  to  evaluate  Pp- 

To  illustrate,  let  the  dimensional  parameter  TV  =  4.  Figure  1  shows  the  performance  of  the 
single-pulse  detector  as  a  function  of  K.  As  K  increases  from  8  to  20,  ROC  approaches  the  known 
covariance  case.  Kelly  has  shown  that  ROC  is  only  a  functiori  of  L.  For  example,  the  case  where 
N  =  16,  K  =  32,  and  L  =  17  yields  the  same  results  as  the  preceding  case  with  TV  =  4  and  K  =  20. 

2.5  Multiple-Pulse  Test  in  a  Stationary  Noise  Enviroiunent  (Case  2) 

The  true  covariance  matrix  is  assumed  to  be  constant  over  the  Np  pulse  observation  intervals 
as  does  the  mean  value  of  the  primary  data  vector.  Analysis  begins  by  forming  an  TV  x  TVp  data 
block  Zt<f  where 


Zn  =  [ziZ7-- -^Nf]  (45) 

The  columns  of  Zs  are  the  primsiry  data  vectors  Zi,  i  =  1, . . . ,  TVp  . 
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Figure  1.  Single*puUe  ie$i  ROC. 


An  Nf  X  Vp  unitary  matrix  U;v,  can  be  appropriately  selected  to  rotate  the  mean  of  the  data 
solely  Into  the  first  column  of  Z//.  This  operation  Is  possible  because  all  the  mean  vectors  are  the 
same.  Performing  the  postmultlpllcatlon  of  Zn  by  results  In 

Zat  =  [Za  Zfc)  ,  (46) 

where  Zo  is  a  column  vector  with  mean  6p  and  Z^  is  an  x  {Np  -  1)  zero-mean  array. 

The  problem  has  now  been  reduced  to  the  case  1  form  with  a  single  primary  vector  with  a 
mean  of  6p  and  K  +  Np  -  I  secondary  vectors.  Test  performance  can  be  evaluated  using  the  Pfa 
and  Pd  equations  develo|jed  In  Section  2.5. 

2.6  Summary  of  Results 

For  the  single-pulse  case,  the  decision  rule  derived  from  a  GLR  test  has  been  statistically 
characterized  under  both  hypotheses.  Under  hypothesis  Ho,  the  PpA  is  solely  a  function  of  the 
excess  number  of  secondaries;  therefore,  the  test  Is  a  CFAR  test.  The  expression  for  the  Pd 
is  complicated  and  requires  considerable  computation  time.  The  ROC  plots  indicate  that  test 
performance  approaches  the  matched  filter  levels  as  the  the  number  of  secondary  data  vectors  used 
in  the  test  Increases. 
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3.  NONCOHERENT  ADAPTIVE  DETECTION 


This  section  examines  the  more  dlfficul  nf’.tatlonary  multiple-observation  case.  Due  to  the 
nonstationarity  assumption,  the  covariance  mav  u  :an  vary  over  the  complete  observation  interval. 
The  true  covariance  matrix  at  each  time  Instant  is  assumed  unknown  under  both  hypotheses  and, 
for  this  case  only,  the  model  for  the  returned  signal  amplitude  is  time  varying.  A  GLR  test 
using  data  from  Np  data  observation  periods  is  then  formulated  under  these  assumptions.  The 
resulting  decision  rule  is  shown  to  be  equivalent  to  a  noncoherent  detection  process.  The  resulting 
test  consists  of  the  Inverse  product  of  Np  j3-distributed  r.v.’s  and  is  independent  of  the  unknown 
covariance  parameters,  implying  that  it  possesses  the  desirable  CFAR  property  [7].  An  exact 
expression  for  the  Pfa  obtained  and  an  upper  bound  on  the  Pp  is  derived  [13).  Theoretical  and 
simulation  results  are  compared  and  observations  with  respect  to  performance  are  made. 

3.1  Multiple-Pulse  GLR  Test 

This  section  formulates  the  GLR  test  for  the  multiple-observation  case,  where  the  covariance 
matrix  and  the  returned  signal  amplitude  are  unknown  parameters.  To  restate,  the  pdf  of  the 
primary  vector  for  the  ith  pulse  under  Hi  is 


(47) 


where 


Mu  =  E  {(z<  - 

bi  is  the  complex  returned  signal  amplitude,  and  p  is  the  steering  vector.  The  primary  vector  pdf 
under  Ho  is 


po(z<|Ho)  = 


(48) 


where 


Moi  =  Eziz/^ 

The  secondary  vectors  are  defined  to  be  zero-mean  and  possess  the  same  pdf  as  that  of  the  primary 
vector  under  Ho. 
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As  per  the  single-puiss  case,  the  joint  pdf  of  the  primarj'  and  secondary  vectors  are  formed 
at  time  instant  i.  Both  likelihood  functions  are  the  product  oi  K  +  l  multivariate  Gaussian  pdf’s. 
Assuming  the  covariance  matrices  are  totally  unknown  under  both  hypotheses,  and  the  returned 
signal  amplitudes  are  unknown  parameters  only  under  Hi,  the  likelihood  ratio  is 


A 


TT  A  —  TT  ^0 

‘”iLi 


(49) 


where 


x<  =  (zi  -  bi  p) 


and 


Qi 


According  to  the  GLR  test  procedure,  the  numerator  of  Equation  (49)  must  be  maximized 
with  respect  to  each  Mh  and  6j.  The  denominator  is  maximized  solely  with  respect  to  each  Moj- 
Taking  the  (/f  +  l)th  root  of  the  numerator  and  denominator  pdf’s  for  each  A*  and  applying  the 
inner  product  identity  x^Y“'x  =  Tr{Y~‘xx^},  the  likelihood  ratio  at  i  is 


where 

S<  =  53 

Ai  is  redefined  as  the  {K  +  l)th  root  of  the  test  at  i. 
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Prom  Equation  (50),  while  fixing  bi  temporarily,  the  MLEs  of  the  covariance  under  each 
hypothesis  are 


Mu  =  ((*i  -  bip)  (zi  -  bip)^  +  Si}  (51) 

and 

M<k  =  7^  {«,"+&}  ■  (52) 

Again,  the  sample  covariance  matrix  estimator  form  arises  from  the  current  problem  formu¬ 
lation,  not  from  an  a  priori  assumption.  After  substituting  the  MLEs  of  the  covariance  matrix  for 
each  pulse  Interval  into  their  respective  likelihood  functions,  the  follovring  ratio  of  determinants, 
which  still  contains  the  unknown  parameter  bi,  is  obtained: 


I|M(h|| 

IIMull 


(53) 


Equation  (53)  is  maximized  with  respect  to  bi  by  minimizing  the  denominator,  i.e.. 


i«l  minllMull 
bi 


(54) 


Using  the  same  determinant  identities  found  in  Section  2.1, 

(A  -H  1)"^ IIMull  =  IlSillll  -h  (z.  -d.p")'^Sr*(2i  -  bip)]  (55) 

and 

(K  -K  l)^|lMo<ll  =  IIS4I11  +  z«Sr*Xi)  ,  (56) 

and  completing  the  square  with  respect  to  the  bi,  an  expression  for  the  GLR,  where  the  unknown 
parameters  6*,  t  =  1, . . . ,  TVp  are  isolated,  is  obtained: 
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Np 

A=n 


1  +  z,^S~4i 


min^  l+Z,^Sj  *Z<+p^S“*p 
bi 


p"sr‘zi 


p^S~*Zi 


(57) 


Equation  (57)  is  maximized  when  bi  =  bi,  i  =  1, . . . ,  Np,  as  was  the  case  in  the  single-pulse  analysis. 
From  the  previous  equation,  the  MLE  for  in  is 


P^S-‘z, 


(58) 


The  estimator  hi  possesses  a  form  equivalent  to  the  prediction  error  (see  Section  2.2).  Substituting 
hi  in  place  of  bi  yields  the  final  form  of  the  LRT: 


A  = 


i+if'S; 

r^zi 

1  +  z,^S"‘zi  - 

p«S-‘zi  ^ 

p'^S-'p 

> 

<  a 


(59) 


where  a  is  a  threshold  parameter.  The  above  expression  is  just  the  product  of  Np  single-pulse  tests. 

The  GLR  test  is  solely  a  function  of  the  primary  vectors,  the  secondary  vectors,  and  the  known 
vector  p.  Taking  the  logarithm  of  the  GLR,  a  summation  of  Np  Independent  Ai  terms  is  obtained. 
Because  the  data  at  each  time  instant  are  used  to  make  decisions,  the  overall  detection  process  is 
equivalent  to  noncoherent  Integration.  An  examination  of  LRT  statistics  given  by  Equation  (59) 
follows. 


3.2  Analysis  of  the  Noncoherent  GLR  liest 

This  section  examines  the  test  statistic  A  under  each  hypothesis.  The  inverse  test  statistic 
under  hypothesis  Ho  possesses  a  central  distribution.  The  cumulative  distribution  function  of 
the  Inverse  test  statistic,  Pfa,  is  evaluated  in  terms  of  a  closed-form  expression.  Under  hypothesis 
Hi,  an  exact  expression  for  the  pdf  of  the  Inverse  product  of  Np  Independent  noncentral  0  r.v.’s 
is  unobtainable.  Hence,  Pp  cannot  be  expressed  in  closed-form  and  an  upper  (Chemoff)  bound  is 
obtained. 

The  equivalence  between  the  problem  studied  here  and  the  well-known  Wilks’  A  test  (10]  is 
discussed  next.  The  general  form  of  the  Wilks’  A  test  arises  when  the  covariance  matrix  is  unknown 
and  the  subsequent  GLR  procedure  results  n  the  ratio  of  determinants.  Through  the  Judicious  use 
of  several  identities  found  in  multivariate  statistics  and  probability  theory,  the  moments  of  the  test 
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can  be  shown  to  be  equal  to  the  moment  of  the  product  of  independent  /3-distributed  r.v.’s  (for 
real  and  complex  r.v.’s,  see  Anderson  (16]  and  Box  [21],  and  Kelly  and  Forsythe  [5],  respectively). 
When  the  product  consists  of  greater  than  three  elements,  however,  an  exact  expression  for  the  pdf 
is  virtually  impossible  to  obtain.  Adding  to  the  complexity  of  this  problem  is  the  presence  of  the 
loss  factor  p,  which  necessitates  working  with  conditional  distributions. 

Many  statistical  researchers  have  examined  the  unconditional  Wilks’  A  test.  Several  asymp¬ 
totic  approaches  based  on  the  distribution  have  been  developed,  but  these  methods  often  require 
using  a  correction  factor  [15,22].  These  approximations  produce  reliable  results  only  when  Np  is 
relatively  small.  Because  this  analysis  may  be  extended  to  cases  where  the  pulse  train  can  be  quite 
long,  the  approaches  do  not  appear  to  be  useful  in  the  current  context.  In  terms  of  notation, 
the  statistical  community  defines  the  likelihood  ratio  pjHo/pjHi,  which  is  the  inverse  of  the  defi¬ 
nition  commonly  used  by  communications  engineers:  pjHi/p]Ho.  The  difference  arises  from  what 
each  group  is  testing.  The  goal  of  the  communications  engineer  is  to  find  differences  between  two 
hypotheses  given  the  data,  whereas  the  goal  of  the  statistician  is  to  determine  whether  the  data  fit 
a  specified  model.  The  test  definition  in  this  report  results  in  the  inverse  product  of  independent 
/?  r.v.’s. 

A  single  term  of  the  GLR  test,  Ai,  given  in  Equation  (59),  was  shown  in  Section  2  to  be 
statistically  equivalent  to  1  plus  the  ratio  of  two  independent  r.v.’s,  \vi\^/Ti.  Under  hypothesis 
Hi,  possesses  a  complex  noncentral  x^  distribution  with  1  complex  dof  and  a  noncentrality 
parameter  Cj  =  Oipi,  and  is  simply  a  complex  central  x*  r.v.  under  H©.  The  noncentrality  parameter 
c^  is  the  product  of  the  nonadaptlve  SNR  Oi  and  a  /?-di6trlbuted  loss  factor  Pi.  Hence,  under  Hi, 
luj]^  is  a  complex  noncentral  x*  r.v.  that  is  conditioned  on  pi.  The  parameter  7<  is  an  unconditional 
complex  central  x*  r.v.  with  L  complex  dof.  In  terras  of  these  x^  parameters,  the  GLR  test  for  the 
problem  under  study  is 


> 

<  a 


(60) 


The  complex  central  x*  pdf  is  defined: 

/x3(x<;  m)  =  , 

where  m  is  the  dof  parameter  and  the  noncentral  x^  pdf  is 

43(x<;m;c<)  Im-i{2y/x^  , 

where  d  =  OiPi  is  the  nonceutrallty  parameter  and  Im-i  is  the  modified  Bessel  function. 


(61) 


(62) 
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Because  the  ratio  of  the  two  Independent  ^  r.v.’s  possesses  an  F  distribution,  the  term  luil^/Tj 
is  conditionally  F-distributed.  Using  the  previously  defined  relationship  between  0  and  F  r.v.’s, 
the  multiple-pulse  GLR  can  be  expressed  as  the  inverse  product  of  ^  r.v.’s  that  are  conditioned  on 
the  lass  factor  pi. 

Under  the  null  hypothesis,  the  general  form  of  the  complex  Wilks’  A  test  can  be  expressed  as 


xp{m  +  1  -  i,  n)  <  a  (63) 

«=i 

The  GLR  test  in  the  noise-only  case  is  a  special  form  of  the  Wilks’  A  test  in  that  the  inverse  of 
each  At  possesses  the  same  distribution,  i.e.,  the  current  test  is  of  the  form 


Np 

n 


1 

xp(m,n) 


> 

<  a 


(64) 


The  conditioning  effect  and  the  constant  dof  parameters  are  cases  not  found  in  the  general  statistical 
literature.  The  distribution  of  the  GLR  under  each  hypothesis  is  now  examined. 

3.2,1  The  Pfa  of  the  GLR  Test 

Pfa  is  defined  P(A)|Ho  dA.  The  numerator  term  Vi  is  assumed  to  be  zero-mean  (therefore, 
the  noncentrality  parameter  =  0),  which  leads  to  the  result  that  the  distribution  of  the  test  is 
a  function  of  central  r.v.’s.  For  the  dof  specified  in  this  problem,  the  GLR  test  statistic  under 
hypothesis  Ho  in  terms  of  the  product  of  central  r.v.’s  is 


1 

X0(L,  1) 


(65) 


Using  the  definition  of  the  complex  /3  pdf,  the  pdf  for  each  Ai  reduces  to 
/^(x;L,l)  = 


(66) 


By  allowing  the  returned  amplitude  to  be  time  varying,  the  single-pulse  case  results  have  been  used 
to  this  point;  however,  the  distribution  of  the  Inverse  product  of  Np  r.v.’s  must  be  determined. 

In  general,  the  distribution  of  the  product  of  r.v.’s  is  extremely  difficult  to  evaluate  analyt¬ 
ically.  If  a  transformation  can  be  performed  so  that  the  test  consists  of  the  sum  of  independent 
r.v.’s,  several  probability  identities  can  then  be  applied  to  obtain  the  desired  pdf  [6].  Letting 
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Ai  =  «p  (I) 


(67) 


and  performing  the  standard  change  of  variable  transformation  proccdv-:e,  the  pdf  of  the  newly 
defined  r.v.  yi  is 


fiVi)  ~  . 


(68) 


Hence,  the  term  yi  is  a  central  r.v.  with  one  complex  dof.  The  test  can  now  be  expressed  as 


Np  /N, 


A  =  Ai  =  exp  I  ^  yi/Lj  =  , 


(69) 


where 


Taking  the  logarithm  of  Equation  (69) ,  the  resulting  test  consists  of  the  sum  of  Independent 
r.v.’s;  therefore,  Y  is  central  x^-dlstrlbuted  with  Np  dof.  The  pdf  of  Y  is  defined: 


/(r)  =  4,(y;iVp)  = 


nNp) 


(70) 


The  cumulative  distribution  function  (cdf)  of  a  x^  r-V-  with  Np  dof  is 


F{Y)  =  l-Gf,,{Y)  , 


(-) 


where 


N,-i 


r4.(i:i/V,)di  =  e-''  £■  ^ 

•'r  ■ 


This  expression  is  the  incomplete  Gamma  function,  and  the  Pfa  of  the  test  can  be  expressed  as 


PfA  ~  Gn JYo) 


(72) 
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where 


yb  =  LlogQ 


The  threshold  that  yields  a  particular  PpA  value  can  be  found  by  solving  Equation  (72)  by  means 
of  a  Newton-Raphson  iteration  procedure.  The  required  derivative  is  given  by 

=  •  (73) 

An  exact  expression  for  the  Pfa.  he-,  the  error  under  hypothesis  Ho,  has  been  obtained. 
The  Pfa,  corresponding  to  a  pulse  train  of  arbitrary  length  is  easily  evaluated  using  the  previous 
expressions.  To  complete  the  analysis,  the  test  statistic  under  Hi  is  now  examined. 

3.2.2  A  Chernoff  Bound  for  the  P[)  of  the  GLR  Test 

Pd  is  defined  /“  p(A)|Hi  dA.  Using  the  same  notation  as  in  the  Section  3.2.1,  the  GLR  test 
is  composed  of  the  inverse  product  of  Np  noncentral  P  r.v.’s.  The  pdf  of  a  single  noncentral  P  r.v. 
with  L  and  1  complex  dof  and  noncentrality  parameter  Ci  is  (5) 


Unfortunately,  no  simple  transformation  can  be  performed  to  enable  the  product  of  terms  to  be 
expressed  in  the  preferable  sum-of-terms  form  as  was  done  in  Section  3.2.1.  As  a  result,  the  problem 
must  be  approached  In  a  different  way. 

Statisticians  have  used  asymptotic  x*  approaches  to  obtain  approximate  solutions;  however, 
many  of  these  approaches  have  drawbacks,  including  the  use  of  correction  factors.  This  study 
chooses  to  obtain  the  Chernoff  bound  for  the  distribution  of  the  LRT  [23).  Although  not  exact,  this 
upper  bound  will  at  least  provide  a  measure  of  confidence  in  results  obtained  through  simulations. 
The  remainder  of  this  section  is  devoted  to  obtaining  an  optimal  (tightest)  upper  bound  on  the 
Pd  for  the  Inverse  product  of  Np  independent,  noncentral  P  r.v.’s.  The  derivation  leading  to  the 
Chernoff  bound  Involves  several  probability  and  Integral  identities  [5,9]. 

Because  the  nonccntrallty  parameter  of  each  r.v.  is  conditioned  on  a  /3-distributed  loss  factor, 
the  unconditional  distribution  of  the  test  statistic  must  be  obtained.  The  GLR  under  Hi,  expressed 
as  the  inverse  product  of  conditional  noncentral  p  r.v.’s,  is 


i-I  ial 


1 

1;  c») 


(75) 
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The  pdf  and  the  moments  of  an  r.v.  that  are  composed  of  the  sum  of  independent  r.v.’s  can 
be  determined  using  moment  generating  function  (MGF)  properties.  To  take  advantage  of  these 
properties,  the  following  change  of  variable  is  made; 

Ai  =  e"'*'  (76) 

The  GLR  test  can  now  be  expressed  as 

N,  Np  /Np  \ 

^  =  II^r‘ =  11®“’ =  (77) 

»=1  i-\  \i=l  / 

Taking  the  logarithm  of  Elquation  (77)  and  defining  U  —  log  A,  the  log  GLR  test  is 

U  =  '^Ui  <\ogQ  .  (78) 

i=l 


Equation  (76)  has  the  form  that  enables  using  MGF  properties. 

Basing  the  remaining  statistical  analysis  on  the  r.v.  U,  an  upper  bound  on  Pr(f/  >  logo  is 
determined.  Setting  6  =  loga,  the  Chemoff  bound  in  terms  of  U  is  defined; 

Pt{U  >6)<  e-^^Ele^^^l  .  (79) 

The  optimization  parameter  u  that  yields  the  tightest  possible  upper  bound  is  the  solution  of 


du 


=  0 


(80) 


To  obtain  the  optimal  upper  bound  for  this  problem,  the  term  E[e*'^),  which  is  the  MGF  of  the 
test  statistic,  and  the  parameter  v  need  to  be  evaluated. 

The  Chemoff  bound  can  obviously  be  expressed  in  terms  of  MGF,  based  on  its  definition.  The 
r.v.  U  is  the  sum  of  independent  r.v.'s;  therefore,  using  the  MGF  summation-product  relationship, 
the  MGF  of  U  is  the  product  of  the  MGFs  of  each  individual  u,; 


isl 


im\ 


(81) 
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Citing  the  definition  of  the  r.v.  iti,  E[e*^]  is  equal  to  E(Ai  *'].  The  MGF  of  Uj  is  clearly  equivalent 
to  the  — i/th  moment  of  a  noncentral  0  r.v. 

The  -J/th  moment  of  the  r.v.  A  is  obtained  by  taking  the  expected  value  of  -i4h  power  of 
Aj.  The  pdf  of  a  noncentral  complex  0  r.v.  with  L  and  1  complex  dof  is  repeated  for  clarity: 


i0{Xi\L,  l;Ci) 


r(L+l) 
r(L  +  1  +  k) 


Ci^'e  ‘^*‘f;j(ii;L,  1  +  k) 


Given  the  pdf  of  Aj,  the  —  t^h  conditional  moment  of  this  r.v.  is 


X(L\  Lr(L-^)  > 

SWr(i--  +  *  +  *)^ 


X  e  ^iF\(l  +  k\L-u  +  l+k;Ci) 


(82) 


(83) 


The  unconditional  moment  of  Aj  is  now  obtained.  Conditioning  is  removed  by  averaging 
over  the  distiibution  of  the  loss  factor.  First  set  d  =  Oipi,  then  multiply  Equation  (83)  by 
fffiPuL  +  l,A’  —  I),  and  finally  integrate  over  the  range  of  pi.  The  resulting  expression  for  the 
— i/th  unconditional  moment  of  Aj  is 


rf/H  =  ^ r(L4N)  TjLik)  , 
/iwr(Llfc-i/)  TiLNk)  r(L  +  l)“’ 

^r(l  +  A:+j)r(LlA:  +  j)  T  (LNk)?  (Llk  ~  u) 

""  ^  T{\  +  k)r{Lik)  r{LNk+j)riLik-u  +  j) 

r- 

where  Llk  =  A  +  1  +  fc  and  LNk  =  L  +  N  +  k. 


(84) 


The  function  iFi{N  ~  1;  L  +  N  -F  k  +  j]ai)  is  the  confluent  hypergeometric  function.  It  will 
converge  to  a  finite  value  because  the  second  parameter  is  greater  than  the  first  (the  denominator 
gamma  function  values  will  grow  faster  than  those  of  the  numerator)  [24].  The  summation  over  the 
index  j,  which  is  of  a  form  similar  to  a  2F2  hypergeometric  series,  will  also  converge  for  the  same 
reason  given  earlier. 

A  compact  expression  for  the  MGF  of  t/  is 


(85) 
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Assuming  that  the  nonadaptive  SNR  is  constant  over  the  observation  interval  (remember  that  only 
the  model  for  the  returned  amplitude  is  time  varying,  not  the  actual  returned  amplitude),  the  MGF 
of  U  reduces  to 

Mt;(i/)  =  (E[Ar])^'’  ■  (86) 

Equations  (84-86)  are  used  for  the  computer-based  evaluation  of  Mt/(i/).  The  final  form  of  the 
expression  for  the  Chemoff  bound  requires  that  the  optimization  parameter  u  be  found  numerically. 
Computational  difficulties  encountered  in  obtaining  the  ChernofT  bound  in  Section  3.3. 

3.3  Simule^ion  Results 

This  section  discusses  a  comparative  analysis  that  is  performed  using  results  obtained  from 
Monte  Carlo  simulations,  the  Chernoff  bound,  and  the  benchmark  matched  filter  case  where  the 
covariance  matrix  is  assumed  known.  How  the  Chernoff  bounds  were  obtained  numerically  is 
described  first,  followed  by  a  discussion  of  the  Monte  Carlo  simulation  technique  used  for  this 
investigation.  Detector  performance  is  then  analyzed  using  ROC  plots.  For  the  purposes  of  this 
report,  ROCs  are  defined  as  plots  of  Pd  versus  SNR  for  a  fixed  PfA'  The  ROC  for  the  known 
covariance  case  is  produced  using  a  modified  Marcum-Q  function  [25]. 

The  numerical  evaluation  of  the  Chernoff  bound  involves  an  expression  that  includes  two 
Infinite  series.  The  \F\{N -\\L+N confluent  hypergeometric  series  is  well  behaved  and 
converges  rapidly.  The  summation  on  i  was  truncated  after  the  difference  between  two  successive 
terms  was  <  10*^.  The  summation  on  j  ,  however,  is  not  quite  that  simple.  This  summation  is 
similar  in  form  to  a  jFj  hypergeometric  series  [24]  that  monotonlcally  increases  to  a  peak  value  and 
then  decreases  until  the  summation  reaches  a  steady-state;  therefore,  truncation  bounds  were  not 
applied  until  after  the  maximum  functional  value  hod  been  reached.  It  was  necessary  to  evaluate 
the  logarithms  of  the  terms  in  the  Chernoff  bound  expression  to  avoid  overflow.  Also,  recursive 
relationships  were  used  to  decrewe  the  computation  time.  For  each  SNR  value,  the  optimization 
parameter  v  was  evaluated  using  a  numerical  search  method  with  v  restricted  to  the  range  0  <  i/  <  L 
due  to  the  definition  of  the  F  function.  The  Chemoff  bound  was  then  computed  for  this  particular 
value  of  v  and  stored  for  plotting  purposes. 

A  simple  Monte  Carlo  approach  was  chosen  to  evaluate  Pp.  For  the  noncoherent  decision 
rule,  it  was  necessary  to  generate  Np  noncentral  ^  r.v.’s,  form  their  Inverse  product,  and  compare 
it  to  a  threshold.  If  the  product  was  greater  than  the  threshold  (corresponding  to  a  target  present 
in  the  data),  a  counter  was  Incremented  by  one. 

A  PpA  of  0.01  was  selected  and  the  corresponding  threshold  was  easily  computed  using  the 
incomplete  gamma  function.  This  relatively  high  (for  radar  applications)  Pfa  value  was  chosen  to 
minimize  the  number  of  simulation  runs.  If  the  decision  rule  is  used  in  conjunction  with  an  M  out 
of  N  detector,  Pfa  '"dll  be  greatly  reduced  [26].  If  one  wanted  to  work  with  Pfa  on  the  order  of 
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10“®,  importance  sampling  or  extreme  value  theory  techniques  could  be  applied  to  the  problem 
under  investigation  (27-30).  Given  the  value  of  the  current  Pfa.  at  least  1,000  data  sequences  of 
Np  length  were  required  for  the  simulation  results  to  be  statistically  significant.  The  final  counter 
value  was  divided  by  1,000  with  the  resulting  value  defined  as  the  simulation-generated  Pd.  This 
process  was  repeated  for  each  SNR. 

The  generation  of  r.v.’s  possessing  the  desired  noncentral  0  distribution  is  now  presented.  A 
noncentral  0  r.v.  is  a  function  of  the  noncentrality  parameter  (c  =  SNR  p)  and  the  dimensional 
parameters  of  the  problem,  L  and  N .  The  nonadaptive  SNR  was  chosen  to  range  between  0  and 
20  dB  in  1-dB  increments  for  simulation  purjjoses.  The  loss  factor  pi  was  set  to  its  mean  value, 
E(Pi)  =  (L  4-  \)/{L  +  N).  Because  pi  is  a  /^-distributed  r.v.  and  its  distribution  is  sharply  peaked 
about  the  mean  for  the  given  dof  parameters,  the  use  of  the  mean  value  in  place  of  the  r.v.  is 
justified. 

For  each  value  of  the  noncentrality  parameter,  the  random  variates  were  generated  using 
a  stratified  rejection  m£*hod  (14,31).  The  rejection  (also  referred  to  as  “acceptance-rejection") 
method  was  selected  because  the  noncentral  0  probability  distribution  is  not  amenable  to  direct 
inversion.  Direct  inversion  requires  solving  X  =  F“*(f/),  where  F(x)  is  the  cumulative  distribution 
function  of  the  r.v.  X,  and  U  is  uniformly  distributed  between  0  and  1.  An  examination  of  the 
noncentral  0  probability  distribution  indicates  that  F~^{U)  would  be  diflScult  to  obtain  for  this 
distribution;  however,  the  rejection  method  can  be  used  with  any  pdf.  The  stratification  process 
significantly  reduces  the  number  of  attempts  required  to  generate  the  desired  quantity  of  r.v ’s.  The 
resulting  computational  savings  is  Important  considering  the  number  of  variates  that  were  required 
for  this  study. 

Analysis  of  detector  performance  capabilities  proceeds  by  examining  the  ROCs  for  various 
parameter  combinations.  Figure  2  illustrates  a  t3T>ical  relationship  between  the  Chernoff  bound, 
the  Marcum-Q  function,  and  the  Monte  Carlo  simulation  results. 

Although  it  is  somewhat  loose  (4  dB  greater  them  the  simulation  results),  the  Chernoff  bound 
provides  some  level  of  confidence  in  the  Monte  Carlo  simulation  results.  In  all  cases  examined, 
the  bound  is  greater  than  the  simulation  result  but  tracks  the  shape  of  both  the  Marcum-Q  and 
simulation  curves.  Notice  that  the  simulation  result  is  only  1  dB  less  than  the  optimal  detector 
over  the  complete  SNR  range.  The  GLR-based  detector  with  N  =  4  and  K  =  20  performs  well 
under  these  conditions. 

The  effect  of  the  excess  number  of  secondaries,  L  =  K  -  N  +  I,  \b  now  examined.  Figure  3 
portrays  the  asymptotic  behavior  of  the  detector  as  a  function  of  K  with  number  of  array  elements 
N  fixed.  As  the  number  of  secondary  vectors  K  Increases,  detector  performance  clearly  approaches 
the  known  covariance  case,  depicted  as  the  Marcum-Q  curve,  with  the  loss  compared  to  the  optimal 
detector  being  significantly  reduced.  The  detector  provides  excellent  performance  given  a  reasonable 
number  of  secondary  data  vectors.  The  Improved  performance  as  function  of  increasing  K  occurs 
for  all  pulse  train  lengths. 
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Figure  2.  Noncoherent  performance  comparison  for  Np  =  2 


Figure  S.  Noncoherent  performance  as  a  function  of  K  for  Np 
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Figures  4  and  5  portray  the  same  information  as  Figure  2  but  with  Np  =  4  and  10,  respectively. 
The  parameter  L  has  the  same  value  {K  =  20,  therefore,  L  =  17)  as  used  for  Figure  1  to  make 
valid  comparisons.  The  ChernofF  bound  tracks  the  Marcum-Q  and  simulation  results  in  both  cases. 
For  both  Np  =  4  and  10,  the  detection  loss  compared  to  the  matched  filter  is  the  same  1-dB  loss 
found  in  Figure  2. 


Figure  4-  Noncoherent  performance  oompariion  for  Np  =  4- 


Figure  6  illustrates  detector  performance  as  a  function  of  Np.  Notice  the  pronounced  leftward 
shift  of  the  ROC  ctirves  as  a  function  of  increasing  Np.  The  SNR  required  to  achieve  a  certain  Pp 
is  significantly  decreased  as  the  pulse  train  is  increased  in  length.  These  results  are  not  unexpected 
because  more  information  is  being  used  to  make  the  decision  with  a  longer  pulse  train. 

To  further  examine  the  effect  of  increasing  Np,  the  SNR  required  to  achieve  a  Pp  =  0.9  versus 
Np  for  =  4  and  =  20  is  shown  in  Figure  7.  The  required  SNR  for  a  fixed  Pp  value  decreases 
at  a  rate  proportional  to  1/ y/T^,  which  is  consistent  with  the  behavior  of  noncoherent  detectors. 
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Figure  7.  SNR  reguired  to  obtain  Pp  =0.9  for  N  =  4  N  =  20. 


It  is  evident  that  the  two  parameters  that  control  detector  performance  are  the  pulse  train 
length  Np  and  the  number,  of  excess  secondary  vectors  L.  The  array  size  N  is  now  increased 
to  evaluate  its  effect.  The  results  illustrated  in  Figure  8  were  produced  under  the  conditions 
N  =  16,  K  =  32,  i.e.,  L  =  17  and  Np  =  2.  The  curves  are  essentially  identical  in  form  to  the  results 
portrayed  in  Figure  2  but  with  L  —  17  and  Np  =  2.  Again,  as  L  was  increased,  the  performance 
approached  that  of  the  matched  filter  in  the  limit. 

The  results  obtained  clearly  follow  the  patterns  established  for  the  single-pulse  case  examined 
in  Kelly  [7].  As  the  number  of  excess  secondaries  is  increased,  a  marked  improvement  occurs 
in  performance.  In  addition,  the  detector  possesses  the  desirable  CFAR  characteristic  in  that  it  is 
statistically  independent  of  the  true  covariance  and  is  only  a  fiinction  of  the  dimensional  parameters 
N  and  K.  Given  that  the  data  can  be  processed  fi'oro  a  large  number  of  secondaries  and  a  long  pulse 
train,  a  detector  that  approaches  the  matched  filter  performance,  even  under  low  SNR  conditions, 
has  been  developed. 

3.4  Summary  of  Results 

A  noncoherent  detector  that  is  statistically  identical  to  a  special  form  of  the  Wilks’  A  test 
has  been  presented.  The  detector,  the  output  of  which  is  the  inverse  product  of  Np  complex  P- 
distributed  r.v.’s,  has  been  analyzed  under  both  the  signal-plus- noise  and  noise-only  hypotheses. 
Under  hypothesis  Ho,  the  detector  is  statistically  equivalent  to  a  central  r-v.  with  Np  dof. 
Under  hypothesis  Hi,  an  exact  expression  for  the  pdf  of  the  LRT  cannot  be  obtained,  and  a 
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Figure  8.  Noncoherent  performance  comparison  for  N  =  16  and  K  =  S2. 


Chernoff  bound  on  the  Pp  was  derived.  The  Chemoff  bound  enabled  some  degree  of  confidence 
in  the  results  obtained  through  Monte  Carlo  simulation  runs.  The  detector  was  shown  to  provide 
excellent  performance  given  that  the  number  of  excess  secondaries  is  greater  than  at  least  bN. 
As  the  number  of  pulses  is  increased,  detector  performamce  improves  as  expected  and  less  SNR  is 
needed  to  achieve  a  certain  Pp  level.  The  use  of  additional  signal-free  data  vectors  was  shown  to 
provide  a  marked  performance  improvement  in  a  detector  derived  from  a  GLR  test  formulation. 
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4.  COHERENT  ADAPTIVE  DETECTION 


An  approximate  solution  foi  the  coherent  detection  problem  is  developed  in  this  section. 
Assuming  a  constant  model  for  the  returned  amplitude,  the  maximization  of  the  likelihood  ratio  is 
mathematically  unmanageable.  As  a  result,  an  approximation  using  a  truncated  series  expansion  of 
the  GLR  forms  the  basis  for  the  coherent  decision  rule;  however,  the  exact  statistics  of  the  coherent 
decision  rule  cannot  be  obtained.  Both  the  threshold  value  corresponding  to  a  specified  PpA  and  the 
ROC  are  evaluated  using  Monte  Carlo  simulations.  Under  most  conditions,  the  coherent  detector 
outperforms  the  noncoherent  known  covariance  case  as  well  as  the  noncoherent  unknown  covariance 
detector  presented  in  Section  3. 

4.1  An  Approximate  Multiple-Pulse  Decision  Rule 

In  Section  3,  the  model  for  the  returned  signal  amplitude  was  allowed  to  vary  as  a  function  of 
time.  Although  it  does  not  match  actual  conditions  over  short  time  intervals,  this  model  enabled 
the  maximization  of  the  GLR  with  respect  to  the  returned  amplitude  to  be  mathematically  man¬ 
ageable.  Noncoherent  detector  performance,  however,  cannot  exceed  the  known  covariance  case, 
the  performance  of  which  is  evaluated  using  the  Marcum-Q  function.  The  goal  is  to  develop  a 
detector  that  uses  the  complete  Np  pulse  observation  Interval  to  obtain  the  MLE  (or  a  reasonable 
approximation)  of  the  returned  amplitude.  An  estimator  of  this  form  enables  a  coherent  decision 
rule  to  be  obtained  '  an  end  result. 

Following  the  ;  '  analysis  used  in  the  formation  of  the  noncoherent  GLR,  the  Np  -pulse 

likelihood  ratio  is  formed: 


pi(zj,Zj(l),...,Zi(K');Mu,fr) 

Po(Zj,  Zi(l), . .  .  ,Zt(R)i  Mot) 


(87) 


where  the  returned  amplitude  6  is  a  constant. 

The  likelihood  functions  in  Equation  (87)  are  first  maximized  with  respect  to  the  unknown 
covariances.  After  taking  the  {K  +  l)th  root  of  the  test  statistic,  substituting  the  MLEs  of  the 
covariance  matrices  into  their  respective  likelihood  functions  and  performing  some  algebraic  ma¬ 
nipulation,  GLR  is 


mlnllMnll 

b 


(88) 


which  is  still  a  function  of  the  unknown  returned  amplitude.  Expanding  the  determinants  in  the 
preceding  GLR  and  conipletlng  the  square  with  respect  to  6,  results  in 
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A=n— 7 

«=i 

mini 

b 

where 

s<  =  f;2i(fc)z«(fc)  . 

isl 

Observing  Equation  (89),  the  maximization  of  the  GLR  with  respect  to  6  requires  the  solution  of 
a  2  Np-degree  polynomial. 

A  suboptlmal  solution  technique  suggested  by  Kelly  (8,12]  begins  by  defining  a  new  parameter 
hi: 


1  4-zf^S~*2i 


i+z,^Si  ‘zj+p"Sj  ‘p 


Hc-K 


)"s-‘ 


Zt 


(89) 


_  P"S-‘z, 
"  P^s- 


In  terms  of  hi,  the  GLR  Is 


_ l-i-zfS:-^z< _ 

inin|i+zfS”'zi+p^^S,~*p|fr-6i(’-p^S7*p(S,|j 


(90) 


(91) 


With  the  definitions 


A  .  l+z,WS-*Zi 

^ - tt: — : — rr 


1 


p  sr‘p 


and 


Vi  = 


P^sr'p 


1  +z^Sf*Z4  ~ 


p^S"‘zj 

7^ 


(92) 


(93) 
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the  likelihood  ratio  can  be  expressed  as 


Np  ,  2^  _l  Ap  Ap  ,  i2\“* 

A=n^oi (' =n^o*n(i+^‘r"^* ) 

i=i  ^  t»i  t=j ' 

At  each  time  Instant  i  throughout  the  observation  lnter\'al.  the  term  hi  is  equivalent  in  form 
to  the  minimum  variance  estimator  of  the  true  returned  amplitude.  Assuming  reasonable  SNRs 
(or  alternatively,  high  quality  covariance  estimates),  the  difference  6  -  hi  should  be  relatively  small. 
Given  this  assumption,  the  GLR  expressed  in  Equation  (94)  is  expanded  using  the  following  identity. 
The  series  expansion  for  a  fraction  of  the  form  (1  +x)~*  is 

-4-  =  1-x  +  x2-...  (95) 

If  X  is  small,  the  series  can  be  truncated  after  the  first-order  term  without  much  loss  in  accuracy: 
therefore,  equating  x  to  b  —  hi,  the  GLR  can  be  approximated  by 

Ap  Ap  .  .  _j  Ap  (  \  i2\ 

i=l  i=l  ^  ^  ,=1  \  »=1  / 

The  approximate  result  in  Equation  (96)  must  still  maximized  with  respect  to  the  unknown 
parameter  b.  Only  the  summation  is  a  function  of  6;  therefore,  concentration  is  solely  on  the  term 
Si^i  ~  ^1^-  Defining  the  term  Soib)  as 

Ap 

SB(h)  =  £v;i<»-6<P  ,  (97) 

i-l 


the  magnitude  squared  term  is  expanded  to  obtain 


(Np'iNp 

SB{b)  =  |6|2 '£Vi-2Re\b*'£ VA  >  +  X; Vi  hi 

i«l  ^  i>=l  J  i»l 


Completing  the  square  with  respect  to  h  yields 


NS  S 

Seib)  =  |h  -  hp  ^  V  +  53  Vi\hi\^  -  |hp  ^  V 

ial  isl  i^l 


(98) 


(99) 
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where 


6  = 


N, 

EK 


The  estimator  6  is  a  normalized,  weighted  sum  that  Is  an  approximate  MLE  of  the  returned 
amplitude.  In  a  stationary  environment,  the  expected  value  of  each  K  would  be  equal;  therefore, 
under  stationary  conditions  the  expected  value  of  the  estimate  would  be  the  average  of  the  Np 
independent  minimum  variance  unbiased  estimates  of  the  returned  amplitude.  The  K  terms  may 
be  viewed  as  weights  that  are  proportional  to  the  SNR.  In  a  nonstationary  environment  the  coherent 
estimator  will  place  a  greater  emphasis  on  estimates  obtained  from  high  SNR  conditions  as  opposed 
to  conditions  where  there  is  much  interference. 

By  substituting  b  for  b  in  the  expression  for  Ssib),  the  minimum  value  is  obtained  and  is 
denoted  as  S^: 


Et’,  K 


As  a  result,  the  final  form  of  the  likelihood  ratio  is 


A  = 


(100) 


(101) 


4.2  The  Statistical  Analysis  of  the  Coherent  GLR  Test 

The  coherent  decision  rule  bears  some  resemblance  to  the  noncoherent  decision  rule.  Each 
Aoi  of  Ekjuation  (101)  is  the  likelihood  ratio  term  found  in  Section  3.  Hence,  in  statistical  terms 
the  coherent  likelihood  ratio  is  the  product  of  (1  —  Sb)  and  the  product  of  Np  /?-distributed  r.v.’s. 
Because  the  coherent  likelihood  ratio  contains  another  level  of  complexity  beyond  the  noncoherent 
likelihood  ratio,  l.e.,  the  presence  of  (1  -  Ss),  it  is  not  expected  that  analysis  will  yield  an  exact 
statistical  characterization  of  the  test.  For  the  sake  of  completeness,  the  statistics  of  Sb  are  now 
examined. 

Following  standard  multivariate  statistical  analysis  procedures  and  steps  presented  earlier,  the 
primary  and  K  secondary  vectors  are  whitened  and  rotated.  Coordinate  rotation  is  accomplished 
via  a  unitary  transformation  that  places  the  nonzero  mean  component  of  the  primary  data  vector 
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into  the  vector  first  element.  The  whitened,  rotated,  Gaussian  N-vectors  defined  in  terms  of  the 
partitioned  components  are 

p  =  [l  0...0]  ,  (102) 

Z<  =  [ZA  Zb]^  , 


Si  = 


S>u  Sab 
Sb/1  Sbb 


and 


Pax  Pab 
Pba  Pbb 


where 


Pi  =  S"‘ 

The  matrix  P»  is  the  inverse  of  K  times  the  sample  covariance  formed  from  the  secondary  vectors 
at  time  instant  t.  The  partitions  of  Zi  ,S~‘ ,  and  Pi  are  Implicitly  understood  to  be  functions  of 
time.  The  subscript  i  is  suppressed  to  make  subsequent  formulas  easier  to  read. 

Next  Sb  Is  expressed  in  terms  of  the  partitions.  (A  more  detailed  description  of  the  terms 
that  follow  can  be  found  in  Section  2.)  For  clarity  the  r.v.  Ti  is  again  defined: 

Tj  =  =  P4A  =  Saa  ~  SabS^^S^a  (103) 

Ti  is  an  unconditional  central  r  v-  with  L  =  K  ~  N  +  \  dof.  The  complex  amplitude  estimator 
hi  possesses  the  form 

hi  =  ZA  ~  SabS^^zb  (104) 

The  parameter  hi  is  itself  a  conditional  Gaussian  r.v.  and  is  the  prediction  error  of  za  given  zb- 
The  /3-distrlbuted  loss  factor  pi  is  defined: 
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(105) 


Pi  = 


1 

1  + *25^^70 


Using  these  definitions,  Vi  in  terms  of  the  partitions  (see  Equation  (95)]  is 

V  =  =  — 

l  +  zgSg^ZB  Ti 


(106) 


which  is  conditionally  distributed  as  the  inverse  of  a  x*  r.v.  Expressing  Sb  (see  Equation  (102)]  in 
terras  of  the  partitioned  components  results  in 


2 

2^  Ti 

[l+ZgS  g^gZa)^ 

isl 

r. 


i=l 


or  equivalently,  as 


isl 


where 


(107) 


(108) 


(109) 


Section  2  showed  that  each  |vi|*/7i  term  is  conditionally  F*distrlbuted;  therefore,  the  first 
summation  in  Sb  in  the  sum  of  Np  conditional  F  r.v.’s.  The  second  term  in  Sb  is  examined  in 
two  parts.  The  denominator,  where  each  Vi  is  equivalent  to  pi/Ti,  is  conditionally  distributed  as 
the  sum  of  the  inverse  of  a  x^  r  v-  The  individual  terms  in  the  numerator  summation,  vSi,  are 
equivalent  to 
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(110) 


(Z^-S^flSp'pZe) 

Vibi  = 

Ti 

Each  term  in  the  summation  is  the  ratio  of  a  conditional  normal  r.v.  and  a  r  v.  The 

overall  numerator  term  is  the  square  of  the  sum  of  Np  such  terms.  Unfortunately,  the  resulting 
distribution  is  not  of  a  common  form-  The  distribution  of  the  ratio  of  the  square  of  this  sum  and 
the  denominator  term  must  also  be  determined.  The  second  part  of  Sb  is  practically  impossible  to 
analyze. 

The  exact  statistical  characterization  of  the  coherent  decision  rule,  let  alone  expressions  for 
the  Pd  and  Ppa,  is  virtually  impossible  to  obtain;  therefore,  the  performance  of  the  decision  rule 
is  evaluated  experimentally,  using  Monte  Carlo  simulations. 

4.3  Simulation  Results 

Because  statistical  analysis  of  the  coherent  decision  rule  is  so  cumbersome,  Monte  Carlo  simu¬ 
lations  were  used  to  obtain  performance  measures.  The  threshold  value  corresponding  to  a  specified 
PpA  was  experimentally  obtained  and  then  used  in  the  evaluation  of  the  Pd-  The  performance  of 
the  coherent  decision  rule  is  analyzed  in  both  stationary  and  nonstationary  environments.  Coher¬ 
ent  detector  performance  (assuming  an  unknown  true  covariance  matrix)  is  shown  to  exceed  the 
performance  of  a  detector  derived  from  the  noncoherent,  known  covariance  matrix. 

The  coherent  decision  rule  upon  which  the  subsequent  results  are  based  is 

Np 

A=  (l-SB)nAo<<Q  ,  (111) 


where  a  is  the  threshold  parameter.  In  Section  4.2,  the  exact  statistics  of  Sb  proved  to  be  unob¬ 
tainable;  therefore,  the  test  statistic  must  be  evaluated  using  the  raw  data  Zi  and  Zi{k)  which  are 
multivariate  Gaussion  random  vectors.  Using  Equations  (4),  (6),  (7),  and  (15),  the  decision  rule 
can  be  computed  simply  in  terms  of  Zi,  Zi(k),  which  is  used  to  evaluate  S~*,  and  the  steering  vector 
P- 

Because  the  steering  vector  enters  the  decision  rule,  a  particular  array  geometry  must  now 
be  assigned.  A  uniformly  (half- wavelength)  spaced  linear  array  was  chosen  for  the  sake  of  compu¬ 
tational  simplicity.  The  N'-element  steering  vector  for  this  array  geometry  is 

p  =  (112) 
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Again,  for  simplicity,  the  target  is  assumed  to  be  located  broadside  to  the  array.  The  steering 
vector  under  this  assumption  has  the  following  form: 

P  =  (l  1...11  .  (113) 

Without  changing  the  test,  p  is  redefined  as  a  unit  vector. 

The  noise  environment  is  composed  of  directional  interference  in  addition  to  Gaussian  back¬ 
ground  noise.  Assuming  that  the  jith  directional  interference  source  is  narrowband,  local  d  in  the 
far-field,  and  that  its  power  is  denoted  as  Oj,  the  total  noise  covariance  matrix  can  be  represented 
as 


Mi  =  12  -  (ll'l) 

j=i 


where  Nj  is  the  number  of  interference  sources,  ffj  is  the  interference  power,  and  is  the  steering 
vector  of  the  interference. 

The  vectors  z,  and  Zi(k)  both  share  the  covariance  matrbe  Mj.  To  create  multivariate  Gaus¬ 
sian  random  vectors  with  covariance  M<,  the  zero-mean  N-vectors  Zu,i  and  Ztn(k)  are  formed  first, 
using  a  random  number  generation  program.  Because  the  elements  of '!  esc  data  vectors  are  sta¬ 
tistically  Independent,  the  resulting  covariance  matrix  is  the  identity  macrix.  1  he  data  vectors  are 
then  statistically  transformed  to  possess  the  desired  covariance  using  the  Choi  sky  factorization  of 
the  true  covariance  matrix.  The  Cholesky  factor  Cj,  which  is  defined  Cj^Ci  =  Mi,  is  applied  to 
the  data  in  the  following  manner: 

Zf  —  CiZ<j,i  (1.15) 

Zi{k)  =  CiZ^iik)  .  (116) 


The  covariance  matrix  of  the  respective  transformed  data  vectors  is  now  the  desired  covariance 
M<.  For  Gaussian  white  noise  processes,  the  transformation  using  the  Cholesky  factorization  is  not 
necessary. 

Under  Hi,  the  expected  value  of  the  primary  vector  Zi  is  bp.  To  form  a  nonzero  mean 
primary,  bp  is  added  to  the  random  primary  vector.  The  value  of  b,  which  corresponds  to  the 
target  radar  cross-sectional  area,  is  computed  assuming  that  the  covariance  matrix  is  the  identity 
matrix.  The  known  white-noise  case  is  used  as  a  reference  for  two  reasons.  First,  assuming  that  the 
noise  is  white  with  unit  variance  and  the  steering  vector  is  of  unit  length,  the  ideal  nonadaptive  SNR 
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reduces  to  |6|*.  Second,  the  reference  SNR  is  now  independent  of  jammer  locations  and  powers.  As 
a  result,  we  can  maintain  constant  SNR  while  varying  the  noise  environment  over  an  observation 
interval;  therefore,  all  subsequent  tests  use  the  same  target  returned  power  for  a  particular  SNR, 
which  results  in  valid  comparisons.  The  mean  returned  amplitude  value  is 


6  =  10^ 


(117) 


where  the  SNR  is  in  decibels.  At :,  the  matrix  S~’  is  obtained  by  forming  the  outer  product  of  K 
secondary  vectors,  summing  the  outer  products  together  and  computing  the  inverse.  The  discussion 
of  the  creation  of  all  the  constituent  components  of  the  decision  rule  is  now  complete. 

As  previously  stated,  the  threshold  cannc.  be  explicitly  calculated  and  must  be  evaluated 
numerically.  In  the  discussion  of  the  simulations  used  to  obtain  the  performance  characteristics  of 
the  noncoherent  detector,  the  reasons  for  the  choice  of  Pfa  =  0.01,  although  very  high  for  radar 
systems,  were  explained.  Assuming  that  the  overall  detection  scheme  incorporates  M  out  of  N 
decision  rules,  the  resulting  Pfa  can  be  brought  down  into  acceptable  levels  for  radar.  Also,  the 
number  of  computations  required  to  produce  statistically  significant  results  is  drastically  reduced 
whei.  compared  with  the  required  number  of  tests  conducted  under  extremely  small  Pfa  conditions. 

The  PfA  is  computed  under  the  assumption  that  the  primary  data  is  zero-mean.  For  a  given 
noise  environment  structure,  10,000  zero-mean  primary  data  vectors  were  used  as  inputs  to  the 
coherent  decision  rule.  At  each  i  during  the  Np  pulse  observation  interval,  a  new  covariance  estimate 
was  formed.  The  decision  rule  outputs  were  then  used  to  create  a  sample  pdf.  The  threshold  value 
corresponding  to  Pfa  =  0.01  was  then  calculated  using  the  definition  of  the  cumulative  distribution 
function.  This  process  was  repeated  several  times,  and  the  average  of  the  threshold  values  was  used 
to  evaluate  the  ROC. 

As  per  the  analysis  of  the  noncoherent  detector,  the  performance  of  the  coherent  detecUjr 
is  analyzed  as  a  function  of  K  and  Np.  In  addition,  the  coherent  detector  is  studied  under  both 
stationary  and  nonstationary  noise  conditions.  For  all  cases  examined  In  the  sequel,  the  Pfa  was 
set  to  0.01.  The  Prj  at.  a  particular  SNR  was  evaluated  using  1,000  trials. 

The  performance  of  the  coherent  detector  was  evaluated  with  N,  the  number  of  array  elements, 
set  equal  to  4.  The  particular  value  of  N  affects  the  performance  of  thr  detector  only  through  the 
parameter  L\  therefore,  the  results  for  N  =  4  can  be  extrapolated  to  r-.'-ized  arrays  by  varying 
K.  The  added  benefit  of  a  small  N  value  is  that  the  number  of  computat  ion;,  i  significantly  reduced 
due  to  smaller  vector  sizes. 

For  the  current  analj  sis  the  known  covariance  case  Is  used  as  the  benchmark  reference.  The 
performance  of  the  noncoherent  detector  under  known  covariance  conditions  Is  obtained  using 
the  Marcum-Q  function.  The  Marcum-Q  designation  in  the  legends  of  Figures  9-11,  14,  and  17 
corresponds  to  this  case.  The  performance  of  coherent  detector  under  these  conditions  was  obtained 
as  follows.  Because  the  covariance  is  assumed  knovm,  the  secondaries  do  not  provide  any  additional 
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information  and  are  dropped  from  the  likelihood  functions  (which  are  the  product  of  Np  identical 
multivariate  Gaussian  distributions,  where  the  true  mean  value  is  unknown  under  hypothesis  Hi). 
After  some  algebraic  manipulation,  the  final  form  of  the  GLR  test  takes  the  form 


A  = 


> 

<  a 


(118) 


After  performing  a  whitening  transformation  and  a  coordinate  rotation,  the  test  can  be  shown 
under  Ho  to  be  equivalent  to  a  central  t-v-  with  one  complex  dof.  The  threshold  value  for  the 
required  Pfa  is  the  logarithm  of  the  Pfa-  Under  Hi,  the  test  is  distributed  as  a  noncentral  r  v. 
with  one  complex  dof  and  a  noncentrality  parameter  c  =  Np  SNR.  The  cdf  of  this  Rician  r.v.  is 
evaluated  using  the  Marcum-Q  function.  The  ROC  corresponding  to  the  coherent  detector  under 
known  covariance  conditions  is  designated  KNOWN  COVARIANCE  in  the  figures. 

Detector  performance  is  now  examined  in  greater  detail.  It  is  assumed  that  the  noise  envi¬ 
ronment  is  known  at  each  i  whether  it  is  stationary  or  nonstationary;  therefore,  as  long  as  the  type 
of  noise  environment  that  is  present  during  the  tth  sample  interval  is  known,  the  interference  is 
always  nulled  (if  present)  due  to  the  S~’  term  found  in  the  decision  rule.  Detector  performance 
is  essentially  independent  of  the  time-varying  nature  of  the  noise  environment  if  the  covariance 
matrix  is  continuously  updated. 

A  typical  case,  A/p  =  4  and  K  =  20 ,  is  shown  in  Figure  9.  The  coherent  detector  clearly 
outperforms  the  noncoherent  detector  examined  in  Section  3.  In  fact,  the  coherent  detector  requires 
about  1  dB  less  signal  power  in  the  low  SNR  region  of  the  ROC  than  the  known  covariance  matched 
filter  represented  by  the  Marcum-Q  curve.  For  this  value  of  K,  coherent  detector  performance  is 
down  only  J  dB  from  optimal  levels.  The  improved  performance  under  the  stated  conditions 
provided  by  the  coherent  detector  was  found  to  occur  for  all  values  of  Np  that  were  examined. 

From  previous  analyses  the  number  of  secondary  data  vectors  has  been  shown  to  play  a  major 
role  in  the  performance  of  the  GLR  test-based  detectors.  Figures  1  and  3  Illustrate  the  asymptotic 
nature  of  performance  of  these  detectors  os  a  function  of  K.  For  the  coherent  detector.  Figure 
10  shows  that  AT  is  a  controlling  factor.  Only  when  K  exceeded  16  did  the  performance  surpass 
that  of  the  noncoherent,  known  covariance  case.  As  K  Increases  from  20  to  30,  the  ROC  begins  to 
appioach  the  theoretical  maximum.  As  the  number  of  pulses  in  the  observation  interval  is  increased, 
the  number  of  secondaries  needed  to  surpass  optimal  noncoherent  detector  performance  increases 
(see  Figure  11). 

Because  of  integration  gain,  as  tbe  number  of  pulses  in  the  observation  interval  is  Increased 
the  coherent  detector  should  provide  increased  the  performance  (less  SNR  required  to  achieve  a 
specified  Pl')-  Figure  12  Illustrates  the  ROC  for  the  reported  standard  array  size  and  K  ~  30 
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Figurt  10.  Coherent  detector  performance  as  a  function  of  K  for  Np=  4. 
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Figure  11.  Coherent  detector  performance  as  a  function  of  K  for  Np=  10. 


as  a  function  of  Np.  These  simulation  results  follow  the  ROC  curves  of  the  matched  filter  as  a 
function  of  the  number  of  samples,  i.e,,  the  curves  move  upward  and  leftward  as  Np  increases.  The 
observation  that  the  integration  gain  effectively  reduces  the  SNR  required  to  achieve  a  specified  Pd 
is  verified  in  Figure  13  for  AT  ts  4  and  K  —  30. 

The  performance  of  the  coherent  detector,  given  a  wide  variety  of  noise  and  interference 
environments,  is  now  examined.  The  presence  of  the  inverse  of  K  times  the  sample  covariance  in 
the  decision  rule  suggests  that  directional  interference  should  be  effectively  nulled  if  K  is  sufiiciently 
large.  The  effect  of  different  interference  powers  and  the  number  of  jammers  on  coherent  detector 
performance  are  considered.  To  perform  this  analysis,  a  stationary  environment  is  assumed  by 
bolding  the  jammer  power  and  the  number  of  jammers  fixed  over  the  observation  interval.  These 
constraints  will  be  relaxed  and  their  effects  examined.  Figure  14  shows  detector  performance  for 
two  different  jammer-to-noise  ratios  (JNRs)  as  virtually  identical.  JNR  is  defined  as  ffjfcr%.  Under 
both  JNR  scenarios  and  white  noise  conditions,  coherent  detector  performance  is  superior  to  that 
of  the  noncoherent  matched  filter.  Figure  15  portrays  the  effect  of  the  number  of  interferers.  The 
detector  produces  essentially  identical  performance  levels  regeudless  of  interference  structure. 

At  this  point  in  the  analysis  it  is  clear  that  as  long  as  is  sufficiently  large  and  the  statistics 
of  the  noise  environment  are  knowm  at  each  sample  Instant,  the  coherent  detector  outperforms  the 
noncoherent  detector  and  the  noncoherent  matched  filter  detector.  It  should  be  noted  that  the 
actual  interference  directions  do  not  affect  performance  unless  they  lie  inside  the  main  lobe  of  the 
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Figure  12.  Coherent  detector  performance  as  a  function  of  Np. 


Figure  IS.  SNR  Required  for  Pq  =  0.9  as  a  function  of  Np. 
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Figure  14.  Coherent  detector  performance  as  a  function  of  JNR. 


Figure  15.  Coherent  detector  performance  as  a  function  of  number  of  jammers. 
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array  pattern.  For  —  4,  the  main  beam  ranges  between  ±30®.  If  a  jammer  lays  inside  the  main 
lobe,  performance  is  seriously  degraded. 

So  far,  the  ability  to  continuously  update  the  covariance  estimate  has  been  a  basic  premise 
of  this  study.  In  actual  practice  this  procedure  involves  collecting  a  large  number  of  data  vectors, 
forming  outer  products,  and  taking  the  inverse  o{  a.n  N  x  N  matrix.  Much  processing  is  required 
and  could  be  very  expensive.  The  issue  of  a  limited  number  of  covariance  estimates  during  an 
observation  interval  is  now  addressed. 

The  stationary  noise  environment  is  examined,  and  two  extreme  cases  are  considered.  Figure 
16  compares  the  performance  of  the  coherent  detector  using  a  single  and  a  continuously  updated 
covariance  estimate.  In  a  stationary  environment  the  additional  covariance  estimates  provide  little, 
if  any,  additional  benefit.  As  in  case  2,  by  rotating  the  data  via  a  unitary  transformation,  the 
problem  can  be  reduced  to  a  single  primary  vector  along  with  K  x  Np  secondary  vectors.  For  a 
fixed  number  of  array  elements,  however,  a  large  increase  in  the  number  of  secondaries  (greater 
than  20  for  N  =  4)  provides  only  small  increases  in  performance;  therefore,  the  results  illustrated 
in  Figure  16  are  consistent  with  existing  results. 


Figure  16.  Coherent  detector  performance  in  a  Jlationary  enwwnmenc 


The  main  goal  of  this  study  was  to  develop  a  detector  that  would  exceed  noncoherent  matched 
filter  performance  levels  in  stationary  and  nonstationary  noise  environments.  It  remains  to  ejtamine 
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the  performance  of  the  coherent  detector  in  nonstationary  environments  with  a  limited  number  of 
covariance  estimates.  The  blinking  jammer  is  an  extreme  case  of  a  nonstationary  noise  environment 
where  the  jammer  randomly  switches  between  ON  and  OFF.  While  OFF,  the  noise  is  simply  white. 
If  the  covariance  estimate  can  be  continuously  updated  regardless  of  the  jammer  state,  the  coherent 
detector  can  achieve  maximum  performance  because  its  structure  enables  it  to  effectively  null  out 
the  jammer  at  each  time  instant.  If  the  number  of  estimates  is  limited  throughout  the  observation 
interval,  the  performance  is  expected  to  degrade  for  the  following  reasons.  Suppose  only  a  single 
covariance  estimate  is  made  at  the  beginning  of  the  data  collection  period.  If  the  jammer  is  ON, 
it  could  be  nulled  out  at  any  time  during  this  period;  however,  if  the  jammer  is  OFF  during  the 
single  sample  of  the  K  secondaries,  the  sample  covariance  matrix  under  these  conditions  is  not  able 
to  provide  nulling  capabilities.  The  result  is  performance  loss. 

Figure  17  provides  a  comprehensive  review  of  these  possible  condition, «  in  a  nonstationary 
environment.  (For  reference  purposes,  the  curve  corresponding  to  continuous  covariance  updating 
is  included.)  Switching  between  ON  and  OFF  was  uniformly  varied  during  the  observation  interval. 
The  single  covariance  white  noise  curve  corresponds  to  the  case  where  the  estimate  contains  only 
white  noise.  Performance  under  this  condition  is  down  5  dB  when  compared  with  the  reference  at 
PD  =  0.09;  however,  when  the  estimate  includes  the  effect  of  the  jammer,  performance  is  equivalent 
to  that  of  the  reference  case.  Because  the  number  of  noise-only  and  jammer-plus-noise  samples  were 
equally  likely,  it  is  not  surprising  that  the  average  performance  curve  lies  roughly  halfway  between 
the  two  previously  mentioned  cases. 
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Figure  17.  Coherent  detector  performance  under  blinking  jammer  conditions. 
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The  conclusion  reached  from  these  results  is  that  if  the  noise  environment  rapidly  fluctuates, 
it  pays  to  frequently  update  the  covariance  estimate;  also  if  a  sufficient  number  of  secondaries  is 
used,  the  coherent  detector  will  exceed  the  noncoherent  matched  filter  levels  even  though  the  true 
covariance  matrix  is  both  unknown  and  time  varying. 

4.4  Summary  of  Results 

A  coherent  detector  that  is  based  upon  an  approximate  GLR  test  has  been  presented.  The 
approximation  was  required  to  enable  the  maximization  of  the  test  with  respect  to  the  constant 
unknown  returned  amplitude  b.  The  detector  was  manipulated  into  a  form  such  that  the  test  is  the 
product  of  the  noncoherent  decision  rule  and  the  term  (1  -  Sb).  Beyond  a  certain  number  of  sec¬ 
ondaries,  the  coherent  detector  performe«i  very  well,  always  surpassing  noncoherent  matched  filter 
and  noncoherent  detector  performance  levels.  Even  under  extreme  conditions  of  nonstationarlty, 
the  coherent  detector  produced  very  good  results,  especially  when  continuous  covariance  estimate 
updating  was  allowed.  In  conclusion,  the  coherent  decision  rule  yields  excellent  adaptive  detection 
performance  in  both  stationary  and  nonstationary  noise  environments. 
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5.  CONCLUSIONS 


Two  new  multiple-pulse  GLR  test-based  detectors  that  provide  excellent  performance  in  non¬ 
stationary  noise  environments  have  been  presented.  Due  to  the  complexity  of  the  analysis,  little 
work  has  been  done  in  the  analysis  of  nonstationary  systems;  however,  because  nonstationary  noise 
environments  are  quite  typical  in  the  radar  context,  the  results  obtained  in  this  study  are  impor¬ 
tant.  One  of  the  strengths  of  this  analysis  is  that  only  precise  analytical  techniques  were  considered 
as  opposed  to  commonly  used  ad  hoc  methods. 

The  analysis  of  the  coherent  and  noncoherent  detectors  are  extensions  of  the  single-pulse  de¬ 
tector  that  was  completely  analyzed  by  Kelly.  The  single-pulse  decision  rule  was  derived  under  the 
assumptions  that  the  true  covariance  matrbi.  and  the  complex  returned  amplitude  were  unknown 
parameters.  To  overcome  covariance  estimation  losses  that  degraded  detection  performance,  addi¬ 
tional  signal-free  data  vectors  were  incorporated  into  the  decision  rule  that  resulted  in  markedly 
improved  performance.  The  decision  rule  was  shown  to  be  statistically  equivalent  to  the  inverse  of 
a  conditional  /3-distrlbuted  r.v. 

The  multiple-pulse  GLR  test  was  also  formulated  using  additional  secondary  vectors.  Also,  the 
true  covariance  matrix  was  allowed  to  vary  over  the  Np  pulse  sampling  intervals.  Unfortunately, 
the  maximization  of  the  likelihood  ratio  with  respect  to  the  returned  amplitude  proved  to  be 
mathematically  cumbersome,  and  to  cope  with  this  problem  two  different  solution  techniques  were 
considered. 

The  first  solution  was  to  allow  the  model  for  the  signal  amplitude  to  be  time  varying  with 
the  resulting  test  statistically  equivalent  to  the  product  of  Np  /^-distributed  r.v.’s.  By  taking 
the  logarithm  of  the  product,  the  test  Involved  the  summation  of  Np  independent  r.v.’s  being 
compared  to  a  threshold,  hence,  the  noncoherent  detector  nomenclature.  Through  the  judicious 
use  of  a  transformation,  the  threshold  value  was  readily  obtained.  The  exact  analysis  of  the  test 
under  Hi  was  shown  to  be  virtually  Impossible,  but  a  Chernoff  bound  was  derived.  As  the  number 
of  secondary  vectors  used  for  covariance  estimation  increased,  simulations  showed  that  noncoherent 
detector  performance  approached  matched  filter  levels. 

The  second  solution  technique  involved  the  use  of  an  algebraic  Identity  that  enabled  the 
test  to  be  expressed  as  a  series  expansion.  Truncating  the  series  after  the  first  term  allowed  the 
test  to  be  maximized  with  respect  to  the  unknown  returned  amplitude.  Because  tix''  romplete 
observation  Interval  (corresponding  to  coherent  processing)  was  used  to  evaluate  the  maximum 
likelihood  estimate  of  b,  the  detector  was  designated  as  coherent.  The  exact  statistics  of  the 
coherent  decision  rule  were  not  obtained  and  Monte  Carlo  simulation  was  used  to  generate  ROCs. 
As  expected,  the  coherent  detector  surpassed  the  performance  of  the  noncoherent  detector  for 
the  same  number  of  secondary  vectors  and  the  coherent  detector  asymptotically  approached  the 
coherent,  known  covariance  detector  performance  levels.  The  coherent  detector  even  proved  capable 
of  handling  the  highly  nonstationary  blinking  jammer. 


In  conclusion,  both  the  coherent  and  noncoherent  structures  provide  a  framework  to  analyze 
detection  performance  in  nonstationary  noise  environments.  The  results  show  that  given  sufficient 
processing  power  and  a  reasonable  number  of  terms  available  for  covariance  estimation,  detectors 
can  be  designed  to  operate  in  rapidly  fluctuating  environments.  In  particular,  the  coherent  detector 
can  provide  almost  optimal  performance  levels  despite  the  necessity  of  estimating  the  covariance 
matrix  and  that  the  true  noise  statistics  are  time  varying. 
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